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Chapter 1 
Introduction 



In spite of many recent developments, in particular in string theory, the problem 
of quantization of General Relativity is still an open one both from mathematical 
and physical point of view. 

In this thesis we discuss some contributions to the lattice approach to Euclidean 
Quantum Gravity, namely to Regge Calculus and Dynamical Triangulations that 
offer the most natural discretization of General Relativity. This discretization con- 
sists in going from a Riemannian manifolds to triangulations of Piecewise-Linear 
manifolds. In both cases the partition function is a sum over the triangulations 
of Piecewise-Linear manifolds. Each triangulation is weighted by a factor equal 
to the exponential of minus the discretized version of the Einstein-Hilbert action 
(Regge-Einstein action). Moreover the diff-invariant continuum measure on the 
Riemannian structures of a manifold M is replaced, in general, by a DeWitt-like 
measure for the edge lengths in Regge Calculus and by a micro-canonical measure 
in the dynamical triangulations. 

The bulk of this thesis is based on two original contributions which are reported 
in chapter three and four, while in chapter two we introduce very briefly the main 
notions about the model of dynamical triangulations. All the notions that we give 
in this chapter are well known and established in the literature and we introduce 
them since they will be used in chapter three. 

Along the main stream of the connection between Euclidean quantum field theory 
and classical statistical mechanics we introduce the notions of micro-canonical, 
canonical, and grand-canonical partition functions. In particular in two dimen- 
sions we mention different analytical ways in which the number of combinato- 
rial inequivalent triangulations was calculated. Recent results on the estimate of 
micro-canonical partition function in three and four dimensions are summarized. 
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This estimate will be extensively used for the classification of the elongated phase 
of dynamical triangulations in four dimensions in chapter three. 
Successively we give the definition of the two point Green function in the context 
of dynamical triangulations. As a consequence of the definition of the Green 
function, we give the definition of the susceptibility function. From its behavior 
near the critical line we define the string susceptibility exponent ystr- 
In chapter three we shall study the elongated phase of dynamical triangulations in 
four dimensions. We begin by introducing Walkup's theorem which characterize 
the triangulations with the topology of the sphere in four dimensions. We review 
the kinematical bounds which are fixed by this theorem. Furthermore by using the 
expression of the estimated canonical partition function, we stress that the average 
curvature is saturated in correspondence to the kinematical bound of Walkup's 
theorem. This means that for values of kz, the inverse of gravitational constant, 
greater than the value k2 for which the Walkup bound is reached, there is in 
the statistical ensemble of equilateral triangulations of the sphere a prevalence 
of particular triangulations called "Staked Spheres". These configurations are 
the only triangulations for which the Walkup bound is realized. They have a 
simple tree-like structure that can be mapped into branched polymers structures. 
Anyway the map is not one to one in the sense that combinatorial inequivalent 
stacked spheres can be mapped into the same branched polymer. We recognize 
that a stacked sphere fits with the model of a network of baby universes which has 
been formulated from the analysis of the results of the Monte Carlo simulations 
in four dimensions. 

We construct two distinct models of branched polymers and we put them in cor- 
respondence with the stacked spheres by the dual map. This analysis shows that 
the string susceptibility of the stacked spheres is less than 1 . 
At the end we analyse a model taken from the theory of random surfaces and 
adapted to the stacked spheres. The aim of this analysis is to show some evidences 
on the analogy between the stacked spheres and self-energy Feynman graphs rela- 
tive to matrix models of two dimensional triangulations. What we learn from this 
analysis and from all previous considerations is that there is a strong evidence 
that the stacked spheres correspond to a mean field phase in which the string 
susceptibility exponents is ys = ^, so that any attempt of performing a continuum 
limit in this phase will give, even if we assume the convergence of the Schwinger 
functions, a Gaussian measure. 

Recent numerical evidence of a first order phase transition of the model of 
dynamical triangulations in four dimensions and the previous strong evidence of 



a trivial elongated phase suggest that more efforts towards new discrete models 
for simplicial quantum gravity might be required. 

Along this line and following the work of various authors we study Regge 
calculus as a local theory of Euclidean Poincare group in the first order formalism. 
The reason for a first order formalism is both theoretical and technical. With 
respect to traditional Regge Calculus the novelty lies both in its formulation as a 
gauge theory and in the first order formalism. The gauge theory approach results 
mainly in the deficit angle being replaced by its sine. The first order formalism 
has others effect of smoothing out some pathological configurations, like "spikes", 
which might prevent the theory from having a smooth continuum limit. These 
configurations are in fact in the region of large deficit angles where the first order 
formalism and the second order formalism are not equivalent on a lattice. 
We first review and improve some definitions of a previous work on this subject. 
In particular it is stressed how a group theoretical formulation of Regge calculus 
allows to write an action on the dual Voronoi' complex of the original simplicial 
complex which is quite similar to a gauge theory and, more precisely, 

looks like the Wilson action for lattice gauge theory. 
We prove that this action does not depend from the orientation of the Voronoi 
plaquette. 

We formulate a first order principle in which we have two sets of independent 
variables: the normals to the n—l -faces and the connection matrices. The normals 
are considered as the analogous of the n-bein in the continuum theory, and the 
connection matrices as the connection one form in General Relativity. The main 
result of this chapter is that we prove in the case of "small deficit angles" that 
Regge calculus is a solution of the first order formalism. This result is not obvious 
if we vary indipendently the two sets of variables above. 

Then we derive the general field equations for the connection matrices and for 
the normals. We use the method of Lagrange multipliers to take in account the 
constraints of the theory. We propose a method for the calculation of the Lagrange 
multipliers by using the one to one correspondence among the normals to the faces 
of the n simplices and the circumcentric coordinates of the vertices. 
A measure for the path integral for this simplicial theory of gravity is introduced 
and it is shown that it is locally invariant under SO(rL). As a last step we propose 
a coupling of this lattice theory of gravity with fermionic matter. This coupling is 
entirely performed by following the general prescription of the continuum theory. 
In other formulation of discrete gravity (Regge calculus and dynamical triangu- 
lations) the coupling with fermionic matter is usually introduced "ad hoc". In 
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this approach the coupling with fermionic matter is given by considering spinorial 
representation of SO(n). 



Chapter 2 

Dynamical Triangulations 



2.1 Introduction 



In this chapter we shall introduce the basic tools for the simplicial approach []l2l 
[T^ []ni to Euclidean Quantum Gravity (see @ for the main articles on the subject 
and also reference f^) via the theory of Dynamical Triangulations |j6^ 
6^ . A precursor of Dynamical Triangulations has been Weingarten [|5l . Romer 



and Zahringer [|2^ proposed for the first time this model as a gauge fixing of 
Regge calculus. In section one we begin by considering the class of equilateral 
triangulations of Picewise-Linear (PL) manifolds (for a review on PL-manifods see 
m HI) ■which are used in dynamical triangulations. Successively we define the 
action for dynamical triangulations as a restriction of the Regge-Einstein action 
III] to the equilateral triangulations of PL-manifolds. The partition function for 
dynamical triangulations is defined over the ensemble of equilateral triangulations 
of PL-manifolds. 

In this framework we address the counting problem of the number of combina- 
torial inequivalent triangulations and we briefly review the different methods of 
enumeration in the two dimensional case. In three and four dimensions some 
recent analytical results are illustrated 



32 33 



In section |2.21| following |6^] we give the definition of the Green functions 
in the framework of dynamical triangulations. The exposition always follows the 
connection with classical statistical mechanics and at the end we give the definition 
of susceptibility function as a direct consequence of the definition of the grand- 
canonical Green function. 

All this chapter has to be considered as an introduction to the main concepts 
of dynamical triangulations which will be used in the study of the elongated phase 
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in 4-D dynamical triangulations. 

Anyway it is important in our opinion to remark that the whole approach of dy- 
namical triangulations must be considered in the general framework of Euclidean 
Lattice field theory as explained in Q chap. 15. It is crucial in this approach to 
determine the Green functions (the Schwinger functions) and to look for a second 
order phase transition in the parameter space. If we have a second order phase 
transition then the quantum theory of gravity may be defined at the critical point 
by looking for the scaling limit of the Green functions. If this limit exists we ask if 
it fits with the Osterwalder-Schrader axioms flsll [Esl. In the affermative case the 



Riemannian Green functions are the Wick rotated version of the Lorentian Green 
Function (the Wightman functions) [^. A first implementation of the above 
ideas of Euclidean lattice field theory to simplicial quantum gravity has been given 
by Rocek and Williams []12| [|13| . 



It can happen (as in the case of Acj)"^) that the continuum limit gives a gaussian 
measure (free theory), then the theory is trivial. We want to stress that the final 
goal of these theories is to find a non trivial continuum limit. 



2.2 The Model of the Dynamical Triangulations 

The standard rule in Dynamical Triangulations (for a review on the recent results 
see |j2ll] and also @) is to consider all the triangulations of PL-manifolds made by 
equilateral simplices of fixed edge lengths, say a. This implies that the geometrical 
structures are even more rigid with respect to Regge Calculus |jl[] (for a recent 
review on Regge calculus see [^). The set of piecewise-linear maps on these 
simplicial complexes depend only by their combinatorial structures. So that two 
triangulations Ta and T'q are equivalent if there is a piecewise linear map cj) between 
them such that it maps one to one the vertices of Ta into the vertices of J' a in 
such a way that [(^[vi),(^[vj)) is an edge of T'a if and only if (vi,Vj] is an edge of 
Tq and so on for every simplex of any dimension. 

The n—2 simplices are called bones B. The dihedral angle []l| of a n simplex ff"^ on 
a bone is cos^^^. If we indicate by q(B) the number of n-dimensional simplices 
which share the bone B, the Regge curvature ^ on the bone B is 

K(B) =2|^27r-q(B)cos-^^^ V(B) (2.1) 

where V(B) is the n — 2-dimensional volume of the bone. By standard formula we 
know that the volume of a n-dimensional equilateral simplex V(cr"^) is 
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V(CT-) = . (2.2) 

The Regge-Einstein action with cosmological term for dynamical triangulations 
Ta without boundary can be seen as a functional of the form 

where the first sum is over all n-dimensional simplices of the triangulation la- 
We define this two bare coupling constants 



^ ^V(a-)+n(n+l)^^^V(B) 
SttG 1 GnG 

. ^ , (-) 

then the action on a triangulation Tq takes the standard form 

S(k^, Ta) = k^Nn - kn-2Nn-2 (2.5) 

where N^t and ^n-i are the numbers, respectively, of the n and n — 2 simplices of 
the triangulation Ta. 

The partition function on the triangulation Ta of a PL-manifolds PL is defined 

as 

Z(PL, kn, kn-2] = Y_ e-^^^^+^^-^^^-' . (2.6) 

Ta 

Notice that all the triangulations of the PL-manifold PL, have the same weight in 
the path-integral (for a discussion about the measure of dynamical triangulations 
see [|2|) . 



The partition function can be rewritten in the following way 

Z{K, K-i] ='^Y. Y. 1 (2.7) 

n Nn Nn_2 #Ta(PL,Nn,Nn-2) 

where means the sum over all piecewise-flat topologies and Y.#TaiVh'Nr, Nn-2) 
is the number of combinatorial inequivalent equilateral triangulations Ta(PL, Nn, N-n-i) 
of a fixed PL-topology, fixed number Nn of n-simplices and fixed number Nn-2 of 
n — 2-simplices. Now we can consider 
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Pa(PL,Nn,Nn-2) = Y_ 1 (2-8) 

#Ta{PL,Nn,Nn-2) 

as the micro-canonical partition function over the ensemble of the equilateral tri- 
angulations of edge length a with fixed PL-topology. The following formula 



(2.9) 



is the canonical partition function, in which 21 1 
Gibbs measure. Finally 



,l<n-2Nn-2 



plays the role of the 



Z(k,,k,_2) = }^}^e-'^'^^-Z(N^,kn-2) (2.10) 

PL Nn 

is the grand-canonical partition function in which e^'^"'^'^ can be considered as the 
equivalent of the chemical potential. Anyway since even in two dimensions the 
sum over the PL-topologies is divergent ( see ref. [|3| for a brief account), many 
times we will restrict to the only topology of the sphere in every dimension. 

Let's start to examine the two dimensional case. The Dehn-Sommerville equa- 
tions (see [p5[ p. 62 and also pO] p. 80) are: 



N2-N1 



f No 
2Ni 



3N2 , 



(2.11) 



where x C^aiVL] ) is the Euler-Poincare characteristic of the PL-manifold VL whose 
triangulation is TaiVL). As is well known in two dimensions, topological manifolds 
and PL-manifolds are equivalent and the topology, for compact, connected and 
orientable two-dimensional manifolds , is completely classified by the genus g of 
the manifold ||3^. From the Dhen-Sommerville equations p.l2| we can express all 
the components of the f- vector (see p. 78) as function of N2 and x (Ta(PL)], 
that is to say 



N7 

= -Y+xiTaiVL)) 



(2.12) 



2 



The asymptotic number of combinatorial inequivalent triangulations of the 
sphere was calculated for the first time by the mathematician Tutte [HOf. The 
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Figure 2.1: Feynman graph for the matrix model relative to the triangulations. 
Each line has two entries one for the row and other for the column of the matrix 
M. In each vertex there are always three double lines 



underlying idea has been to map by the stereographic projection the triangulations 
of the sphere in planar triangulations and to enumerate them by the techniques 
of generating series built up by using the geometrical properties of the planar 
triangulations. The result is the following 
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Pa(S^N2]xN2^e'^2N2 ^ (2.13) 

where is a numerical constant. It is important that the growth of the number 
of the triangulations is at most exponential. In fact in the opposite case the 
divergences make it impossible to define a statistical theory like ^77| . The result 
|2.13| can be obtained again by using the quantum field theory techniques of the 



matrix models . This techniques is based on the use of the generating functional 



ZN(g) = jdM exp (^-jtr{M^) - -^tr{M')^ (2.14) 

in which M is an N x N Hermitian matrix and A the coupling constant, dM is 
the Haar-measure on this matrix group. The path-integral |2.14| will generate a 
perturbative series whose Feynman diagrams can be represented as double line 
one for each index of the matrix Mij (see figure 2.1). 



Since in the equation |2.14| there is a trace of the product of matrices, the 
Feynman diagrams have to be closed graphs. The fact that there are two lines 
and the possibility of all combinations of the indices means that in general the 
graphs can be closed only on Riemann surfaces (see figure pT2|) . Furthermore each 
Feynman graph corresponds to a dual triangulation. 

It can be show that the perturbative series in A relative to connected Feynman 
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Figure 2.2: An example of a second order perturbation term which is rappresented 
by two double graphs of Tr(M^) which can be closed on a two dimensional torus 



double graphs generated by p.l4| can be arranged as power serie in N of the genus 



g of the Riemann surfaces. It follows that the coefficients of this series give the 
generating functions for the number of combinatorial inequivalent dual triangula- 
tions of a fixed topology. Since there is a one to one correspondence between any 
triangulation and its (topological or metrical) dual we can obtain the generating 
functions of the combinatorial inequivalent triangulations for any topology (see 
111 and also |4|]). 

Another derivation of the formula p.l3| has been recently found and is based 



on the enumeration of the curvature assignment of the bones of the triangulation. 
This estimate can be 

extended to all two-dimensional topologies and to higher dimensions. 

It follows that the two dimensional grandcanonical partition function can be 
written as 



Z(k2,ko] = ^e^^'9'_^Nj^*^"^e-f'^^-^-'^^)N2 (2.15) 
g N2 

where x(g]=2 — 2g. k2 — ^ — k^is called the critical line since for ki > ^ + k^ 
the partition function is convergent and for k2 < ^ + ^ the partition function is 
divergent. jstT is the exponent of the string susceptibility. This exponent, that in 
two dimensions is equal to —j, controls the divergent part of the partition function. 

More precisely if we fix the genus g and consider the limit ki ^ + 1<2) , we 
have that from |2.1E 



lim Zg(k2, ko) = Zl'^ + (k2 - ^ - k^)'-^-^ (2.16) 
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with the finite part (if any) of the partition function near the critical line. 
It is now clear that the discontinuity points of the partition function are on this 
line, so an eventual critical point is on this line. 

As regard the sum over the genus g of the surface it is divergent and even not 
Borel summable. There is a (non-rigorous) way to define a renormalizzation for 
this series known as double scaling limit (see for a brief account). Anyway, in 
general, we will restrict to a given topology, for example the sphere S^. 

In three and four dimensions in the case of the sphere topology the estimate 
|30[] gives, for the microcanonical partition function, an exponential bound. In 



particular by posing N„^_2 = Nnil we have that 30 



Pa(S^N,,N,_2) = g(Nn,Ti)e^-^(^\ n = 3,4 , (2.17) 

where g(NrL,ri) has a subleading growing in respect to the exponential factor 
gNnq(ri) ^j^j ^Yie following geometrical arguments fix the range of r] in 

the closed interval [rii.riJ. So that using the Euler-Maclaurin summation formula 



44] we can approximate the sum by the integral so that the canonical partition 



function is 



■r|2 



Til 



g(Nn,Ti)e^-''i'^''^-2ldTi (2.18) 



where c|i(ri,kTT^_2] = qIti) +f]^n-2- We can use the Laplace method for giving 
an asymptotic estimate of the integral |2.18| . We have to compute the point of 
absolute maximum ri* of the function qi(ri,krL-2) in the interval of integration. In 
general this point will be a function of kn-2, that is r|* = ri*(kn-2)- So we have 

Z(S^N-,k^_2) X g(N^,Ti*(k,_2))e^'^^^'^**'^-^''^^' . (2.19) 

We define k^(kn-2) = qi (■n*(l<^n-2), k2], so that the grand-canonical partition func- 
tion is 

Z(S",k^,kn-2) x^g(Nn,Ti*(k,_2))e-^'^'^"-^"^-^» . (2.20) 

It is clear that now the critical line is k^ = k^{kn-2) with the same meaning as in 
the two dimensional case. 
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2.3 Green Functions in Dynamical Triangulations 

Since the geometry in General Relativity has a dynamical role, the definition of 
the Green function in Quantum Gravity is different from ordinary quantum field 
theory. 

The formal continuum definition of the unnormalized two point function is the 
following 



where A and G are respectively the cosmological and gravitational constant, the 
sum is over the topological structures of the differentiable manifold M, the measure 
V(g) is over the Riemannian structure allowed by M, Se_h is the Einstein-Hilbert 
action and dg{x,y) is the geodesic distance between the points x and y on the 
Riemannian manifold (M, g). To normalize GilA, G,r) we have to divide it by 
the partition function. So in this definition we have to perform a sort of average 
over all points of the Riemannian manifold (M, g) that are at distance r. In the 
discrete we cannot transfer verbatim the above definition due to the coordinate 
invariance of the dynamical triangulations. It is necessary to define the two point 
function in such a way that the analogous quantities, in the discrete, of x and y 
in the continuum have a dependence from the cut-off a. 

A path in a triangulation Ta of a n-dimensional PL-manifolds is a sequence 
of I n-dimensional simplices with the property that Sj and Sj+i (Sj+i ^ Sj) 
have a common face, that is to say I(Sj,Sj+i) = 1 (see p. 526). Suppose that 
the path does not intersect itself, that is to say I(Sp,Sq) = 1 if and only if or 
q=p — 1 orq=p + l,we call it a simple path. In this case we can define the 
length of the path {Sj}-^^ as the number of the faces that the I simplices have in 
common, that are I — 1 . This definition comes from graph theory, in fact if we 
think to the dual of the path {Sjjj^^, every simplex is mapped into a vertex, and 
the face between two simplices is mapped into a edge joining the two vertices. So 
the length of the path is equal to the number of the dual edge of the path (see fig. 
[23|). 




We define the distance between two simplices A and B of the triangulation 
as the length of the simple path which has the minimal length among all simple 
paths {Sj}-^^ such that Si = A and Si^ = B. 



G2(A,G,r)= Y_ V[g)e 



SE-H(M,g] 



Top(M) 




(2.21) 
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Figure 2.3: A simplicial path between the two simplices Si and Sn of a trian- 
gulation. The simplicial path is simple and its dual is obtained by joining the 
baricenters or the circumcenters of the path's simplices. On this path the proper- 
ties of the incidence matrix I are illustrated 
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Let's now give the definition of unnormalized micro-canonical Green function. 



The unnormalized micro-canonical Green function is defined as 64 



#Ta(PL,r,Nn,Nn-2) 



(2.22) 



where H^jafPLrNn Nn-2) number of inequivalent equilateral triangulations 

Tq('PL, r, Nn, NrL-2) of a PL-manifold PL with two labelled n-simplices at fixed 
distance r and with fixed number Nn, and N^-i of respectively n and n — 2 sim- 
plices. The normalized partition function is obtained by dividing p.22| by the 
microcanonical partition function. So on we have that the canonical Green func- 
tion is 



G(PL,r,N,,kn_2)= Y_ e-'^-^^-^G(PL,r,N,,N,_2) 



(2.23) 



and finally the grand-canonical Green function 



G(PL, r, kn, kn_2) = Y_ e''"''"G(PL, r, N,, k,_2) 



(2.24) 



Nn 



Prom the grand-canonical Green Function we can define another quantity often 
used in the literature the susceptibility xi'P^^ ^n, ^n-z) which is defined in the 
following way 



xiVL, kn, K-z) = Y_ ^' 



(2.25) 



r=0 



In the previous definition we can exchange the sum over r with all the sums in 



the definition of the grand-canonical Green function p.24| up the micro-canonical 
Green function |2.22| . So that we will have 



x(PL,k^,k^ 



Nn 



knNn 



Nn-2 



,-kn-2Nn-: 



T #TaCPL,r,Nn,Nn-2) 



1 



(2.26) 



The last two sums in the left hand side of the previous equation are the number of 
inequivalent triangulations with two labelled n-dimensional simplices. Por large 
value of we will expect that asymmetrical triangulations will prevail in the 
micro-canonical ensemble (this is, for example, true for two-dimensional convex 
polyhedra [^, for same particular planar triangulations [Hf and planar maps |^ 



and it is conjectured for all planar maps. In four dimensions this is verified, so to 
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say, indirectly because once we do this antsaz and derive the formula for the num- 
ber of baby universes{see appendix |BD it is in good agreement with the numerical 
simulations (HI). In this hypothesis for each labelling of a n-simplex we 



will produce a number of combinatorially inequivalent triangulation that asymp- 
totically is NnPal'PL, Nn, Nn-2)- Finally the leading asymptotic estimate of the 
number of combinatorially inequivalent triangulations with two labelled simplices 
is 

Y_ Y. xN^pJPL,N^,N,_2) (2.27) 

r #Ta(PL,r,Nn,N,x-2) 

Then the susceptibility p.25| is equal to 

X{VL,K,K_2) X Y_ N^e-'^"^'^Z(PL,N,,k,_2) (2.28) 

Nn-2 

SO that 

X(PL,kn,kn-2) X ^^Z(kn,kn_2,^I-) . (2.29) 

Near the critical line kn i-^ (k^(kTv-2])^ we will have the following asymptotic 
behaviour of the susceptibility 

lim xlVL, kn, kn-2) X x""^ + ^4 (2.30) 

kn^(k5^(kn-2))+^ > n, n ZJ A (kn-k^(kn_2))^=*^ 

where x^^^ is the finite part (if any) of the susceptibility. 
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Chapter 3 

The Geometry of the Elongated 
Phase of 4-D Simplicial Quantum 
Gravity 



3.1 Introduction 



In this chapter we are going to discuss in detail the elongated phase of 4-D dy- 
namical triangulations. The early numerical simulations in four dimensions have 
produced some evidence of the presence of two distinct phases in the ensemble of 



four dimensional triangulations: the crumpled phase and the elongated phase [49 



5^ [0 (it is interesting to remark that this picture emerges also in the simulations 



of 4D quantum Regge calculus []T5| []T6[ [jT^])- The crumpled phase is characterized 
by the presence of few bones and many four simplices incident on them. There 
are very few baby universes (see appendix §). The elongated phase is character- 
ized by the presence of many bones and of many baby universes of smallest size 
(blips). These facts have encouraged people to use the relative abundance of baby 
universes as an order parameter to distinguish these two phases. 
Anyway it was still unclear why in the elongated phase there is an upper kine- 
matical bound on the values of r] = ^ and a prevalence of simple tree-graphs 
(branched polymers) which appear as a model of proliferating baby universes [|9| . 
Established in in section ^]2| we explain that the upper kinematical bound on 
the values of r| is due to the Walkup's theorem for the triangulations of the four 
sphere. In particular this theorem fixes the upper value of r| in the large-N4 limit 
N4 I— > oo. In section |3]3| we introduce a recent [^Of analytical estimate for the 
canonical partition functions for dynamical triangulations in four dimensions rel- 
ative to the triangulations of the sphere S"*. This estimate has the right behaviour 
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in the limit N4 i-^ co, because the average curvature |3.16| in the elongated phase 
remains constant in perfect agreement with Walkup's theorem. Furthermore this 
theorem says also that the dominating configurations in this phase are " Stacked 
Spheres" . In section |3]4| we discuss how the stacked spheres have a simple tree 
structure and can be mapped into branched polymers structures (see appendix 
0). However the problem arises that the map is not a one to one, since there 
are more stacked spheres than branched polymers. So in section ^ we consider 
the action of dynamical triangulations restricted to stacked spheres and show that 
the partition function over the ensemble of the stacked spheres is equal or greater 
than the partition function of branched polymers which are the immages of the 
above map, and equal or less than a particular model of branched polymers. This 
argument shows that the susceptibility exponent ys is less than one. To show 
that y s = ^ we use an argument taken from the physical literature ||6^ . This 
argument follows the line of a model introduced for the first time in the theory of 



random surfaces [163 



We take into account some particular kind of non standard triangulations and 
put it in correspondence with the model of stacked spheres considered as a model 
of proliferating baby universes. Physical arguments ensure that the two models 
belong to the same universality class and we show that the susceptibility exponent 
is ys = ^ as for the branched polymers. 

These types of analytical results were already well known from Monte Carlo simu- 
lations. An early analytical attempt, in the direction of interpreting the elongated 
phase as a stacked spheres phase, has been done in reference |j5^ in which a map 
from a sort of stacked sphere triangulations to branched polymers is considered. 
Anyway the deep reason for polymerization mechanism and with the kinematical 
bound on the values of r\ is given only by Walkup theorem. 

All these facts point out that the elongated phase of 4-D dynamical triangu- 
lations is a mean field phase, so that we expect that if we take the continuum 
limit in the framework of Euclidean lattice field theory (c.f. [|3^ chap. 15), in this 
phase, we will obtain a probabilistic Feynman measure which is Gaussian, or, in 
other words, a trivial theory 



3.2 Walkup's Theorem 

As anticipated in section of chapter || due to geometrical constraints the vari- 
able ri in four dimensions can vary in a finite interval [rii.rji]. The Hmits r[^ and 
r\2 of this interval are fixed respectively by the following Walkup theorem and 
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Dehn-Sommerville equations [|6^ . 

Let us label the bones by an index oc and denote by q(cx] the number of 
four simplexes incident on it; the average number of simplexes incident on a bone 
is 

-^L..^.-0(^) (3.) 

since each simplex is incident on 10 different bones, b ranges between two kine- 
matical bounds which follow from Dehn-Sommerville and a theorem by Walkup 



60|] . This latter theorem is also relevant in classifying the "elongated phase" of 



four-dimensional simplicial quantum gravity as we shall see below 

Theorem 3.2.1 If T is a triangulation of a closed, connected four- dimensional 
manifold then 

Nim>5No(T)-yx(T) (3.2) 

Moreover, equality holds if and only iflE H'*(1 — ^x(T)), where the class of 
triangulations H'^{n) is defined inductively according to: (a) The boundary 
complex of any abstract five- simplex (Bda) is a member o/H'*(0). (b) //K is 
in H'*(0) and a is a four-simplex of K, then K' is in H'^(O), where K' is any 
complex obtained from K by deleting a and adding the join of the boundary 
complex Bdff and a new vertex distinct from the vertices o/ K. (c) If K is 
in H'^(rL), then K' is in H'^(rL+ 1) if there exist two four- simplexes cri and Gz 
with no common vertices and a dimension preserving simplicial map cf) from 
K — CTi — 0^2 onto K' which identifies Bdo"i with Bdcri but otherwise is one to 
one. 



In other words H'^(O) is built up by gluing together five-dimensional simplexes 
through their four dimensional faces and considering only the boundary of this 
resulting complex. H'^(rL) differs from \\^[0] by the fact that it has n handles. 
This way of constructing a triangulation of a four-sphere has a natural connection 
with the definition of a baby universe (see appendix 0). 

A baby universe is associated with a triangulation in which we can distinguish two 
pieces. A piece that contains the majority of the simplices of the triangulation 
that is called the "mother", and a small part called the "baby". In the "minbus" 
(minimun neck baby universes) the two parts are glued together along the bound- 
ary of a four dimensional simplex (in four dimensions) that is the " neck" of the 
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baby universes. Thus the "stacked spheres" can be considered as a network of 
minbus in which the mother is disappeared and the baby universes have a min- 
imal volume, that is the boundary of a five simplices minus the simplices of the 
necks through which they are glued to the others . We will exploit this parallel in 
section to give an estimate of the number of distinct stacked spheres. 
In four-dimensions Dehn-Sommerville equations fen] are 



No-Ni +N2-N3 + N4 = X (3.3) 
2Ni -3N2 + 4N3-5N4 = 
2N3 = 5N4 , 



substituting the third equation in the second we have 



Ni=^N2-^N4 , (3.4) 



and substituting these last equation in the first equation of 3.4 we get 



This imply that 



No-^N2 + N4=X • (3.5) 



IOy 

b<5 + ^ . (3.6) 



and substituting and |3^ in [S^ we obtain 



b>4-15^ . (3.7) 
Thus in the Hmit of large N2, 4 < b < 5. 

3.3 Canonical Partition Function 

Now we introduce more extensively the results in four dimensions for the canonical 
partition function which has been rapidly summarized in section |2]2] of chapter ||. 
Recently it has been analytically shown Q that the dynamical triangulations 



in four dimensions are characterized by two phases: a strong coupling phase in 
the region — log9/8 < k2 < kf^^ and a weak coupling phase for k2 > k2"^^. 
k^'*^* is the value of k2 for which in the infinite volume limit the theory has a 
phase transition from the strong to the weak phase (for a detailed analysis see 
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pod). The transition between these two phases is characterized by the fact that 
the sub dominant asymptotic of the number of distinct triangulations passes from 
an exponential to polynomial behaviour (c.f.||3^). Presently the precise value of 
k2^^* has not been established yet, it is just known that it is close but distinct from 
^max _ i^og4, recent numerical simulations suggest that kj^'^^ = 1 .24. 
In the strong coupling phase the leading term of the asymptotic expansion of the 
canonical partition function is (c.f. pC|] ) 



-4 



Z(N4,k2) = 04 ( '-^^^i^^r:^ ] N4"^exp[-m(Ti*)N-H] 



3A(k2 



exp 



where for notational convenience we have set 

"27 ' 



A(k2 



+ 



2 

27 



,k2 ^ 1 



,27 



e^2 + ]]2-] 



1/3 



1/3 



and 



Tl*(k2 



1 



3^^ ~ A(k2) 



(3.8) 



(3.9) 



(3.10) 



The explicit form of m(ri*(k2)) and tih, the Hausdorff dimension, are at present 
unknown. 

In the weak phase k2 > k^"^* the number of distinct dynamical triangulations 
with equal curvature assignment have a power law behaviour in N4. The phase 
is characterized by two distinct asymptotic regimes of the leading term of the 
canonical partition function: the critical coupling regime and the weak coupling 
regime. In the critical coupling regime k2"* < k2 < k^'^" we have 



Z(N4,k2] X 04 



exp 



^ A(k2)+2 ^ 
. 3A(k2) ) 

[lOlog^ 



-4 



-5 



(3.11) 
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in which T(ri*) is not explicitly known. 

Instead in the weak coupling regime ki > kj^ 



271 /(I _3t^^^J(1 -2ti, 



g(N4+l)f(rimax,k2) 

kf P - k2 

in which rj"^'^'' = ^ for the sphere S^.N^^ 1 = 10(N4 + 1 ] and 



(3.12) 



f(Tl,k2) =-Tll0gTl + (l -2Tl)l0g(l -2ti)-(1 -3Tl)l0g(l -3Tl)+k2Tl . (3.13) 

Prom this form of the canonical partition function it follows that in the case of the 
sphere S"^ and in the infinite volume limit the average value of b, see equation (|3.1|) , 
is a decreasing function of k2 in the critical coupling regime and it is constantly- 
equal to 4 in the weak coupling regime, that is to say 

lim <b>N4=^, kf^<k2<kr'' 

lim <b>N,= 4, k2>k2"'''' . (3.14) 

N4— >oo 

If we look at the average curvature we have 



lim r!- < y K(B)vol4(B) >= ttq^Vs ( 10Ti*(k2) - -cos^^l ) 
N4^oo N4 \ n A J 



(3.15) 



and 



,max 
^2 



1 /5 5 1 \ 

lim — <y K(B)vol4(B) >= na^VZ cos"^- k2 > k^ 

N4^oo N4 B \2 7r 4/ 

(3.16) 

in which vol4(B) is the volume of the four simplexes incident on the bone B, 
and it is a'^^qlB). As we have already said the value of the average curvature 
is saturated. 



This result was already known in numerical simulations [|9[ and, as we will see 
in detail in the next section, it is due to the prevalence of stacked spheres in the 
weak coupling regime. 
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3.4 Elongated Phase 

We have seen that in the case of triangulations of the sphere S"^ for h-) log 4, 
in the infinite volume limit < b >n4*-^ 4. Walkup's theorem Q implies that the 
minimum value of b is reached on triangulations K of the sphere that belong 
to H'^iO) (stacked spheres). Then for ki > log 4 one has < b >n4= 4, that means 
that in this region of kz the statistical ensemble of quantum gravity is strongly- 
dominated by stacked spheres. 

As expHcitly shown in the elements of H'^(O) can be put in correspondence 
with a tree structure. Let us recall that a d— dimensional simplicial complex T, 
d > 1, is called a simple d— tree if it is the closure of its d— simplexes 0i, cTt 
and these d— simplexes can be ordered in such a way that: 

CI o-j fl I U CI c^ij = CI Tj (3.17) 

for some (d— 1 ]-face Tj of ffj, j > 2, and where the Tj are all distinct. This ordering 
of the simplexes of T induces a natural ordering of its vertices in vi, ...,Vt+d, where 
V|_|_d is the vertex of CTi not in CI Ti. Note that the interior part of T contains the 
simplexes CTi and faces Ti. The boundary of T, Bd T, consists of the boundary of 
the fft minus the Ti, and is topologically equivalent to S'^^^ 
It can be shown [pO|, and it is very easy to check, that any element of K G H'^(O) 



is the boundary of a simple (d + l]-tree T, and that K determines uniquely the 
simple (d + 1 )-tree T for d > 2. 

Note that any stacked sphere K G H'^(0] can be mapped into a tree graph (see 
also the reference [^). This mapping is defined in the following way, let's consider 



the (unique) simple five-simple tree T associated with K, every five-simplex is 
mapped into a vertex and every four-dimensional face in common with two five 
simplexes is mapped into an edge which has endpoints at the two vertices which 
represent the two five simplexes (see fig. 3.1). 

Since the map between K and T is one to one, we have a map from a stacked- 
sphere into a tree-graph whose number of links at every vertex can be at most 6 
( since a five simplex has six faces). This map, from the stacked spheres H'^[0) 
to the simple tree-graphs, is not one to one. The mathematical reason is that the 
previous construction maps every T to a tree-graph by an application that is the 
dual map restricted to the five and four dimensional simplexes of T in its domain 
and whose image is the tree graphs that are the 1 -dimensional skeleton of the dual 
complex. It is well known that the dual map is a one to one correspondence only 



28 The Geometry of the Elongated Phase of 4-D Simplicial Quantum Gravity 




if we take in account the simplexes of T of any dimension. It follows that our map 
is such that a simple tree-graph may correspond to many stacked spheres. To 
illustrate better this point we shall use a picture closer to the physical intuition. 
Let us consider a stacked sphere K and its simple tree T. Consider on T one 
of its ordered simplexes CTj and let Tj be the face that it shares with one of the 
simplexes CTi,...,CTt introduced previously. Repeating the same arguments used 
in the calculations of the inequivalent triangulations for baby universes (see 
also appendix 0), we can cut T in Tj in two parts, such that each part has two 
copies of Tj as a part of its boundary. Since Tj is a four dimensional simplex it 
is easy to admit that we can glue this two different copies in 5 -4 -3 ways to 
rebuild again a simple five tree T'. In general all the possible gluings will generate 
distinct triangulations and consequently distinct stacked spheres (if the two parts 
are highly symmetric triangulations some of the 60 ways of joining will not be 
distinct but, since for large TM4 asymmetrical triangulations will dominate, the 
number of case in which this will happen will be negligible). The corresponding 
tree graph associated with these configurations T' will always be the same since 
this operation has not modified the one skeleton of the tree T. 

In statistical mechanics simple-tree structures correspond to "branched-polymers" 
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(see appendix So we have that for k2 > kj^"" the dominant configurations for 
the triangulations of S"^ are branched polymers. This fact was already observed in 
numerical simulations in four dimensions. In [^a network of baby universes of 
minimum neck (minbu) and of minimum size (blips) without a mother universe 
was obtained. These have been interpreted as branched polymer like structures. 



3.5 Stacked Spheres and Branched Polymers 

In this section we will establish a parallel between the standard mean field theory 
of branched polymers and the stacked spheres, in the sense that we will use the 
counting techniques of branched polymers to give an upper and lower estimate of 
the number of inequivalent stacked spheres. As we have stressed in the previous 
section there exists a map between the stacked spheres and tree graphs and this 
map is not one to one in the sense that there are more stacked spheres than tree 
graphs. This means that the number of inequivalent stacked spheres with a fixed 
number N5 of five-simplexes is bounded below by the corresponding number of 
tree graphs. Now we will study the statistical behaviour of the tree graphs using 
the measure of simplicial quantum gravity and restricting it to the stacked spheres. 
In this analysis we shall follow the theory of branched polymers as explained in 
appendix 0. The reader may refer to it for the details. A more mathematical 
analysis on the enumeration of inequivalent tree graphs is contained in [|6^ . 



First of all we notice that the boundary of a five-simplex has six four-simplexes 
and every edge of a tree graph correspond to a cancellation of two four-simplexes 
in the corresponding boundary of the stacked sphere. Since in a tree graph with 
N5 vertices there are IM5 — 1 edges, we have that the boundary of a stacked sphere, 
whose corresponding tree graph has N5 points, is made by a number N4 of four- 
simplexes 

N4-4N5 + 2 . (3.18) 



Prom the condition for stacked spheres (|3.2|) , we have that the Einstein-Hilbert 



action, for dynamical triangulations, restricted to the stacked spheres is 

S = N4(k4- ^k2) -5k2 . (3.19) 

It follows that the Gibbs factor for the ensemble of tree graphs (branched 
polymers) corresponding to stacked spheres is 
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exp (-(4N5 + 2)(k4--k2)+5k2 



(3.20) 



Following the same notation as in [A.l| , let r^(N5) be the number of inequivalent 
rooted tree graphs with N5 vertices and with at most six incident edges on each 
vertex, with one rooted vertex with one incident edge. £,^(N5) is the number of 



inequivalent tree graphs with N5 vertices. By |A.3| we have 



1^5 



(3.21) 



Now let R^(k4,k2) be the partition function for rooted tree graphs with statistical 
weight ( |3.20|) , and Z^(k4,k2) be the partition function for unrooted tree graphs. 
We have 



R^(k4,k2] = Y_ e-(4(N5-l)+2)(k4-|k2)+5k2^6(isj3) 



Ms =2 

00 



Z^(k4,k2) = ^ e-(4N5+2)(l^-|l<,)+51<,^6(^^^ 



Let's define 



(3.22) 
(3.23) 



R*^(Ak4) = (e-^<'^-^'R^(k4,k2)) = e-4(N5-i)Ak4^6^^^^ ^3 24) 



N5=2 

00 



Z*^(Ak4) = (e-^"^-^'Z^(k4,k2)) = e-4N5Ak4£^6(^^ 



N5=1 



where Ak4 = k4 — |k2. It is easy to see, from |3.21| , that 



R*^(Ak4) = -^^^Z*^(Ak4 



(3.25) 



(3.26) 



4dAk4 

Prom |A.11| we know that the asymptotic behavior of the susceptibility x(l<^4> ^2) is 
given by 



Z(k4,k2) 



(3.27) 



Then 
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x(Ak4 



d(Ak4j 



r^(Ak4) 



(3.28) 



V d(Ak4) 
Now we shall study the critical behaviour of this system in order to obtain 



some information about the critical behaviour of the stacked spheres. 
Prom the identity |A.12| we have that 



R*^(Ak4) = 



(|X_l[r^(Ak4)]-j 



(3.29) 



This last equation suggest to define a function F(R*'', Ak4) of R*^ and Ak4 and to 
apply the implicit function theorem to it in 

order to obtain information on Ak4 as a function of R*^. 
More precisely the following function F(R*^, Ak4) of R*^ and Ak4 



F(R*^ Ak4) = R*^ - e"'*^'^ Y_ 



Ky=o 



at the point ((R**^)", (Ak4)°) where 

F((r^]°,(Ak4]°] =0 , 
if it is true the following condition 

aF 



aAk4 



({R*6)0,{Ak4)0) 



defines by the implicit function theorem, locally, Ak4 as 
function of R*^. 



Differentiating equation |3.29| we get a differential equation for R*^(Ak4 



(3.30) 



(3.31) 



(3.32) 



R*^(Ak4) = -4R*''(Ak4 



dAk. 



So R*^(Ak4) shows a singularity when 



g-4Ak4 

1 + (R*^(Ak4))^ - R*^(Ak4] 



H -1 



R*^((Ak4)^) = 1 + 



,-4(Ak4) 



5! 



R*^((Ak4)^; 



(3.33) 



(3.34) 



This equation has only one positive solution and since ■^^^R*^|{ai<4)'= diverges, the 
inverse function tends to zero 
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d(Ak4)(R 



dR 



*6 



^ 



(3.35) 



;Ak4) 



Observe that F(R*^, Ak4) is an analytic function, and in the point ((R*^)", (Ak4)'') 
fixed by equation |3.34| , we have 



3F 



9Akd 



((R*6)S{A1<4)<^) 



\y=o 



(3.36) 



and pairwise 

a? 



aR*6 



= 1 -e 



-4(Ak4) 



({R*6)S(A1<4)<^) 



(3.37) 



this last equation explains why we cannot express in the neighbourhood of (^(R*^)'^, ( Ak4) 
R*^ as function of Ak4 by applying the implicit function theorem to |3.30| . These ar- 
guments are enough to say that in the neighbourhood of ((R**')'^, (Ak4)'^^ (Ak4)(R* 
is an analytic function of R*^. Let us fix the following notation 



)*6i 



Frc = 



a? 



aR 



*6 



^Ak-^ = 



a? 



({R*6)MAk4)'=) 

Frc rc 



a^F 



aAk4 



({R*6)MAk4)c) 



a(R*6)2 

and so on. It is straightforward to calculate that 



({R*6)c,(Ak4)'^) 



p _ p-4{Ak4) 

iR'^R'^ — — e 



)~ ^[(R*')^M <0 



Applying again the implicit function theorem we obtain 



d2(Ak4)(R* 



d(R 



*6]2 



F 



R'^R'^ 



>0 



(Ak4) 



"Ak"; 



All these facts allow to write the following expansion 



(3.38) 



(3.39) 



(3.40) 



Ak4(R*^) - Ak4' 



2 d^(Ak4)(R*^ 
2 d(R*6)2 

+0 ('(R*^(Ak4] - R*^(Ak4 



R*6({Ak4) 



R*^(Ak4) - R*^(Ak4^; 



(3.41) 
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So near (Ak4]'^ we can write 

r^(Ak4) X r^((Ak4)") + C ((Ak4) - (Ak4)")^ (3.42) 
From (|]2|) we get 

x'(Ak4) X ((Ak4] - (Ak4)^)-^ (3.43) 

This means that the critical exponent of the susceptibility for the system of these 
tree graphs is y = ^. 

The partition function that has been studied is a lower estimate of the partition 
function of the stacked spheres. Consider, now, a stacked sphere K and a face Tj 
through which two five-simplexes are glued together (a link on the corresponding 
tree graph). We can glue two four face of a stacked sphere in 5-4-3 different 
ways, for large number of simplexes this will generate distinct configurations of 
stacked spheres whose associated tree-graph is always the same. Repeating the 
same argument for every j, j = 1,...,N5 — 1, (that is to say for every link of 
the dual tree graph) we obtain a factor, (5-4-3)'^^^\ that multiplied by r^(N5) 
and £,''(N5), gives an upper bound on the number of , respectively, rooted and 
unrooted inequivalent stacked spheres. In other words 

Ztree(l<4, ^l] < Z^A^^, ^l) < ltree{^4, ^l) (3.44) 

where Ztree(T<4, ^2] is the partition function for tree graphs studied above, Zs.s.(k4, Vj] 
is the partition function for stacked spheres and ZtTee(l<4, ^z) is the partition func- 
tion for tree graphs with the additional weight defined above. 
A similar analysis as for ZtTee(l<4) ^2) shows that the critical line of Ztreel^A, ^2] is, 
of course, a straight line parallel and above respect to the k4 axis to Ztree(T<4, T<2) 
one, and the susceptibility exponent is again y — \- 



Obviously the estimates (|3.44|) are true for the canonical partition functions too. 



Ztree(N4, ^^2) < Zs.s.(N4, ^^2] < Ztree(N4, ^^2] , (3.45) 

and the previous calculations show that 

Ztree(N4,k2) X N4 ^ e-^^* '^"^^^-^^ ' (3.46) 

and 

Ztree(N4,k2) X N7e-^^f'^-t'^2+t4-llog60) ^ (3^47) 
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where ±4 is a constant that may be calculated by the equation |3.34| . More easily 
from the table in reference we get that ±4 ^ \log 0.34. 

The circumstance that in the weak coupling region the partition function of quan- 
tum gravity, as found in an analytically way in |3^, is strongly dominated by 
stacked spheres, allows us to write 



-3gN4k5{k2) 



(3.48) 



where ys is the susceptibility exponent 



Thus the critical line of the stacked spheres is a straight line parallel and among 
the critical lines of the systems of the two branched polymers. This implies 



k4 



5 5 1 

-k2 + 14 < k5(k2) < k4 - -k2 + 14 - -log 60 



(3.49) 



Moreover from the equations |3.46| and |3.47| the one loop green functions ||64 



of the two model of branched polymers have respectively , near their critical lines 
, the asymptotic behaviour 



Gtree(Ak4) X COSti + (Ak4 - A\C4'']'^ Gtree(Ak4) X COSti + (Ak4 - Ak4'')i 

(3.50) 

This last equation together the equations |3.45| and |3.48| prove that the one loop 



function of the stacked spheres Gs.s.(k4, kz) near the critical line has the asymptotic 
behaviour 

Gs.s.(l<4,k2] xcost3+(k4-k5(k2))^-^^ , (3.51) 

with Ys < 1 (the value ys = 1 is not allowed because in this case the one loop 
green function of stacked spheres at the critical line would have a behaviour like 
^N4=5 ^ ^^^^ divergent and then incompatible with the upper bound given 
by the second equation of |3.50| ). 



Motivated by physical considerations, we can use a well known argument |]64 



57| in favor of the fact that the susceptibility exponent of the stacked spheres is 
ys = ^. More precisely we will show that a model of proliferating baby universes, 
with the measure of quantum gravity restricted to stacked spheres, can be put in 
correspondence with the statistical system of stacked spheres. 
Let us consider four dimensional triangulations that are (boundary of the) stacked 
spheres in which there can be loops made by two three-dimensional simplexes. 
This is possible whenever the stacked spheres are pinched on a three simplex 
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creating a bottle neck loop of two three-simplexes. These loops could be either 
the loops of a Green function or the minimum bottle neck of a baby universes. This 
class of triangulations, following the notation in literature, is called T2. The other 
class of triangulations is the stacked spheres in which the minimum loop length 
can be made by the boundary of a four-simplex that are five three-simplexes. We 
call this last class T5. 

In the two dimensional theory, the introduction in T2 of two-link loops (the 
two dimensional analogous of two three-simplex loop) corresponds in the matrix 
model (|)^ to consider Feynman diagrams with self-energy (c.f. JH^ [|67| ). 
Since the minima loops of T2 and T5, for which they differ, are of the order of 
lattice spacing we will expect that the two classes of triangulations, as a statistical 
mechanics system, coincide in the scaling limit, that is to say they belong to the 
same universality class. 

Let's consider the minimum neck one loop function G(Ak4) in T2. In every 
triangulation of T2 we can cut out the maximal size baby universe of minimum 
neck and close the two three-simplex loop. We will obtain again a triangulation 
that belongs to T2. In this way we will obtain all the triangulations of T2 from the 
triangulations of the stacked spheres T5 considering that for each three-simplex 
either leave them in their actual form or we can open the triangulation to create 
a two three-simplex loop and gluing on it a whole one loop universe G(Ak4). We 
note that in the triangulations of T5 N3 = 5/2N4 (Dehn-Sommerville). CalHng 
the one loop function of T5 G(Ak4), the above considerations lead to the identity 



G(Ak4) = Y_ + G(Ak4])^^ = Y_ e'^'^^ = G(Ak4) , (3.52) 

TeTs TeTs 

in which Ak4 = k4 — 5/2k2 comes out from restricting the action of quantum 
gravity to stacked spheres |3.19| , and where we have defined 



Ak4 = Ak4 - -log (1 + G(Ak4)) 



(3.53) 



By |3.52| we can also write last equation as 



Ak4 = Ak^ + -log (l + G( Ak^ 



(3.54) 



By universality and the estimates (|3.45| ) it follows that near the critical point we 
have that G(Ak4) x cost + (Ak4 - Ak^)^-^^ with Ys < 1 • 
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The susceptibility functions of T2 and T5 by |3.52| are 



x(Ak4 



Prom 3.54 we have 



dAk. 



d(Ak4 



■G(Ak4) x(Ak4) 



dAk. 



:G(Ak4) 



(3.55) 



1 



x(Ak4 



d(Ak4) ■ 2(1+G(Ak4)' ' ^^'^^"^ 

If we calculate di derivative with respect to Ak4 of the one loop function G(Ak4) 
and use the previous equation we get 



x(Ak4 



x(Ak4 



1 



5 x(Ak4) 



(3.57) 



2 (l+G(Ak4) 

Now it is clear that x(Ak4] 1-^ +00 for Ak4 1-^ (Ak4)'^ and with the same 
critical exponent x(Ak4) 1— > +00 for Ak4 h-> (Ak4)'^. Prom equation (|3.57| ) when 
x(Ak4) ^ +00 we have that x(Ak^) ^ §(l + G (Ak^(Ak4'=))) < +00, then the 
system T5 is above his critical line, i.e. Ak4(Ak4'^) > Ak4'^. These facts imply 
that at Ak^(Ak4^) d(Ak4)/d(Ak^) = and around it x(Ak^)/(1 + G(Ak^)) is a 
decreasing monotonic function by equation |3.56| because x(Ak4) 1-^ +co , we can 



expand equation p.54| and |3.57| around Ak4( Ak4 



Ak4 



x(Ak4 



Ak4' 



cost (Ak4- Ak4(Ak4^; 



1 



1 



Ak4^ 



Ak4-Ak4(Ak4^) 

The second asymptotic equality of the last equation implies 



(3.58) 



(3.59) 



1 

2 



(3.60) 



The dominance of stacked spheres in the weak phase allows us to fix the param- 
eter T in the partition function of quantum gravity in the weak coupling regime. 



(3.61) 



T 2 2 

These considerations put evidence that the branched polymer phase of 4-D 
Dynamical Triangulations is a mean field phase. 
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4.1 Introduction 



The model of dynamical triangulations, that we have analysed up to now, has been 
very popular in recent years. The reason of this is due to the excellent results in 
two dimensions where the critical exponents coincide perfectly with the critical 
exponents of two dimensional continuum quantum gravity. This has encouraged 
people to explore 

the model in four dimensions. The main aim has been to explore the dynamics 
of this discrete theory in order to define a four-dimensional quantum theory of 
gravity at the critical point (if any) of the parameters space. At the beginnings 
a second order phase transition has been observed by numerical simulations [|2| . 



Successive and more accurate simulations ||5^ have shown that the phase 
transition is of first order. 

At the moment the situation is not completely clear and further investigations 
seem necessary, but the general agreement of the scientific community about the 
first order nature of the phase transition points out that it is necessary to find a 
new lattice theory of gravity which fits with the general requirements of Euclidean 
lattice field theory. The following model is an attempt in this direction and might 
be an answer to this request. It has been proposed originally and independently by 
different authors ^ |j7|] [|7^, and the original part of the content of this chapter, 
based on the work [^, is a continuation and an improvement of these attempts, 
in particular of [|72| (connected to this work are [^)- We begin in section 
by reviewing some concepts contained in [|7^ and in particular we show how to 
associate to each dual Voronoi edge (see appendix |C| and references lH] and 
p. 395) of the original n-dimensional simplicial complex a Poincare transfor- 
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mation by fixing an orthonormal reference frame in each n-dimensional simplex. 
Then we recall the group theoretical action which is a functional of the sine of 
the deficit angles. Respect to the original work it is shown that the action is 
independent from the orientation of the hinges and from the starting n-simplex 
in which it is written . 

In section ^ following ||7^ we write the group theoretical action 
completely defined on the dual Voronoi complex of the original simplicial com- 
plex. Two sets of variables play an important role: the connection matrices 
U^(a, (3) and the normals b^pla) to the n— 1 faces of each n-simplex a. In the sec- 
ond order formalism the connection matrices Ug(cx, (3) are functions of the b^ptoc). 
As anticipated in we follow a first order formalism in which we consider the 
Ug((x, |3) and the \)^p,[(x) as independent variables. A first order formahsm for 
Regge calculus has been done by J. Barrett []77[ . 



We impose that the first order field variables satisfy only the two constraints 4.16 



and |4.17| . Finally, respect to , in the first order formalism we define a modified 
action in such a way that it is dependent only from each plaquette of the dual 
Voronoi complex . 



In section g]J it is proved that the modified action is independent from the 
orientation of each plaquette. 

In section ^ the field equations for the connection matrices Ug(a, |3) are 
derived in the first order formalism in the approximation of " small deficit angles" . 
It is proved that Regge calculus is a solution for these equations. This result is 
no longer obvious in the first order formalism when we consider the connection 
matrices independent from the b'^^[(x). 

The equivalence of the first order formalism to the (second order) Regge calculus 
is the main new result in this chapter. 

In section we obtain the general field equations on lattice. These equations 
are divided in two sets. The set of equations for the connection matrices U^(a, |3) 
obtained by requiring that the action is stationary under the variation of the 
Ug(a, |3), and the set of equations for b^|3(a). These two sets of equations are 
obtained by using the method of the Lagrange multipliers through which we take 
into account the two sets of constraints, one set on the linear dependence among 
the b^pfa) of the same simplex cx and the second set on the condition that the 
b^l3(a) and bp^((3), as the expressions of the same normal to the n — 1 face seen, 
respectively, in the orthonormal reference frame of a and (3, are connected by the 
matrix U^(a, |3). At the end of section |0| we introduce a method to calculate the 
Lagrange multipliers too. 



4.2 Group Action for Simplicial Gravity 



39 



In section ^ we define a path-integral for this classical theory of gravity on 
lattice and we prove that this quantum measure is invariant under local transfor- 
mation of SO(n). 



Finally in section g]^ we show as it is possible to define in this framework 
a coupling with fermionic matter following, as usual in lattice field theory, the 
example of the continuum theory. 

The main feature of this model is that it is more general than Regge calculus 
and dynamical triangulations. The introduction of a reference frame in each n- 
dimensional simplex is the key concept for which this version of discrete gravity 
is different from the previous one. This is a way to translate on lattice the general 
principle of local invariance under a general Lie group which implies a structure 
, in the continuum, of a Principal Fiber Bundle with the fibers isomorphic to a 
Lie group (gauge theories) or to orthonormal reference frames connected among 
them by matrices of SO(n) (General Relativity in the Riemmanian version). In 
this scheme we can naturally write on lattice the coupling of gravity with other 
fundamental forces and consider extensions to more general theories than gravity 
like, for example, supergravity (see for example WB and [?]). 



4.2 Group Action for Simplicial Gravity 

In this section we shall introduce a group action for simplicial gravity on the 
dual Voronoi complex (see appendix |C| for a definition of Voronoi cells ) of the 
original simplicial complex. In particular we will see that it is possible to associate 
to each dual edge [o^,o^+^) a "Poincare" transformation U(a, a + 1 ) as in lattice 
gauge theories. Of course by a "Poincare" transformation from now on we mean 
in general a SO(n) rotation and a translation. 

We are now going to introduce a reference frame in each n-dimensional simplex 
and we shall see how these reference frames can be connected by using coordinates 
of the vertices of the common n— 1 -dimensional faces and a notion of Levi-Civita 
connection on simplicial manifolds. Consider an hinge h and let {Pi, Pn-i} its 
vertices. Suppose that this hinge is shared by N n-simplices {Si,...,S]m} whose 
vertices are labelled in this way 

S„ = {?!,... ,Pn-l,Qa-l,a,Qa.a+l} (a=1,...,N) . (4.1) 

In each simplex we can fix an origin and a reference frame. In this frame 
the vertices of the simplex have the following coordinates 
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Pi = {yl'la)} a = 1, ...,rL , i = 1, ...,rL- 1 

Qoc-1,a = 

Qa,a+1 = (4.2) 

Let's indicate by the circumcenter of the simplex a and let x'^ioc) their 
coordinates. 

Now it is possible to associate uniquely to a dual edge an element of the Poincare 
group U(a, a + 1] = {U^(cx, a+ 1],U'^(a, a + 1)} by requiring that 



U{:(cc,a+1)y}'((x+1)+U'^(cx,a+1) = y^[oi] 
U^(a,a+1)z^_„+i(a+1)+m(cx,a+1) = z^,,+i(a) , (4.3) 

in other words we are embedding oc and cx + 1 in M"^ and adopting the standard 
notion of parallel displacement in and we transport the origin of the reference 
frame of oc + 1 from a + 1 to the origin of the orthonormal reference frame of a, 
[H |7^. This operation makes the position vectors in cx + 1 of the vertices of the 



common face SanSa+i coincident with the position vectors of the same vertices in 
a. It follows that the matrix U{J(a, a + 1 ) is an orthogonal matrix which describes 
the change from the reference frame of o^+^ to a, considered now as two different 
reference frame of the same vector space ([@]). Since the simplicial manifold is 
assumed orientable [@] we can choose the reference frame in a and a + 1 in such 
a way that Ug(a, a + 1 ) are 
SO(n) matrices. 

The family of all these matrices determine a unique connection that is the Levi 
Civita or Regge connection. 

The arbitrariness of the choice of the reference frame in each simplex oc im- 
plies a sort of local gauge invariance under Poincare transformations A(a) = 
{Ag((x),A-(a)}: 

U(a, a+1) ^ A(a)U(a, a+1)A"^(a + l) , (4.4) 

in particular the coordinates of the dual vertices (the circumcenters of the simplices 
Set) transform as x^{oc) t-> A^{ix)x^{ix) + A°'{ix). Anyway from now on we decide 
to put the origin of the reference frames in the circumcenters, in such a way we 
always have x^[(x) — 0. 
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To the hinge h, it corresponds the dual two-dimensional plaquette that we still 
label by H. Now consider the following plaquette variable 

U|^^J = U(a,a+ l)U(a+ 1,a + 2)...U(a- . (4.5) 

It is evident that leaves the coordinates of the vertices of the hinge h, un- 
changed and its translational part is zero. So it is a rotation of an angle 9(h.) in 



the two-dimensional plain orthogonal to the hinge. It is easy to see that the 
angle Q{h) is the deficit angle [|1| 



N 



e(h) =27r-^e(a,H) , (4.6) 



cc=l 



where 9(a, H) are the dihedral angle of the simplices with the hinge h. 

Let's choose n — 2 linear independent edge vectors — yf — yn-^, •••^^-2 = 
~ Un-1 belonging to h, , and consider the oriented volume V*^' '^^{a] of h, 

V"^' ^^cx) J^^<■■c.-^^^' (aj-.E^-la) . (4.7) 
At this point it seems natural to propose the following gravitational action 

where U^J^^ are the elements of the orthogonal matrix If we choose a 
reference frame in which the first two axis, 1 and 2, are on the two-dimensional 
plane orthogonal to the hinge h the matrix U^y^^ will be diagonal, with the 
diagonal elements equal to 1 , except in the intersection of the first two rows and 
columns where it is 



cos 9(h,) —sin 9(h,) 
sin 9(h,) cos Q{h) 



(4.9) 



Pairwise in this reference frame the only non-zero components of the antisymmet- 
ric tensor V*'^' "^{oc) are V'^' ^^{oc) and V'^' ^^{(X-), because in this reference frame 
the n — 2 independent vectors of h, have zero components on the first two axes. In 
particular the component V*^' ^^((x) of the antisymmetric tensor is equal to V(h,) 
by the definition of the volume V(H) of the hinge H. Since the trace in |0| is 
independent from the orthogonal reference frame chosen in the n-simplex a, the 
action is equal to 
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I = ^sine(h) V(h) , (4.10) 

h 

that for smah deficit angles 9(H) reduces to Regge action 

Notice that in the definition ^ we have chosen the sequence of simplices a, a + 
1 , a — 1 , (X. This sequence induces a direction of rotation in the two-dimensional 
plane orthogonal to the hinge h. determined by the rotation matrix We have 
implicitly assumed that this direction is the same direction generated by the or- 
dered sequence of axes {1 , 2} in the rigid rotation from the axis 1 to the axis 2. 
Since the simplicial manifolds under consideration are assumed orientable @ , we 
can choose the same orientation in each n-dimensional simplex 0. This orienta- 
tion induces an orientation on each face but no orientation on the hinge because 
the two faces of the simplex, that share the hinge, induce opposite orientations 
on it. It follows that in principle there is no preferred orientation on the two- 
dimensional plane orthogonal to the hinge. So we can pairwise choose the se- 



quence a, (X — 1 , CX + 1 , a, the rotation matrix will be the transposed of |475| and 
the direction of rotation is now from the axis 2 to the axis 1 , that is to say the 
opposite of the previous one. In other words respect to the first orientation, now 
the axis 1 ' coincides with 2 and 2' with 1 . It follows that this change of orientation 
implies that e'^'' ^'^ -^'^ = -e^"^^ -^- so that V'^' ^'^V) = V"^' ^^a) . 



Bearing in mind that the analogous of ^ with this new orientation is equal to 



its transposed , we have that the action now ^ is equal again to |4.10| . In this way 
we have proved that the action ^ is independent from the chosen orientation. 
The action ^ by construction is invariant under local Poincare transformations 
(gauge invariance) and is independent from the starting simplex a. This last 
remark may be done more clear by the following lemma 



Lemma 4.2.1 The action ({4-W does not depend from the reference frame 
where it is written, in the sense that if cx and 5t are two simplices which have 
in common the same hinge h then 

(^V'^' a.aAoc] = *^'V(^' a^aA^i) (4-11) 

Dim : Let {a...5i} the simplicial path from the simplex cx to 5t such that the 
simplices of the path share the hinge h. We have that from the definition of 
Levi-Civita connection 



y,-((x)-y^i(cx) = (U((x,a+1)...U(6t_i,50)^(y?(50-y^i(50) (4.12) 
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So for SO(rL] connection matrices is true the following 



V 



(H) 



a, Q2 



[oc] = (U(a,a+1)...U(6t_i,50)::j 



.U(a+1,a))::^ 



(4.13) 



By substituting this last equation in 



(4.14) 



we get the equality (|4.11|) . 

It is interesting to note that this action is similar to Wilson action |j8|] for lattice 
gauge theory (for basic notions of lattice field theory see [|2^ chap. 6). Anyway 
in the form |4]8| the action is in a second order formalism since the plaquette 
variables W^^'^^{h) are compHcated functions by ^ of the coordinates of the 
vertices. Moreover the field variables belong both to the original triangulation 
and its Voronoi dual. We are now going to introduce a formalism that is written 
completely on the Voronoi complex and is a first order formalism like the Palatini 
version [S5| of general relativity. 



4.3 Action on the Dual Voronoi Complex 

Let's consider the n—l -dimensional face f^p = {Pi , Prt} between the n-dimensional 
simplices a and (3. We label the coordinates of the vertices of the face f ocp by (a) 
in the reference frame oc. We define the vector b^g 



= eab,..b,_,E?'...E^V • (4.15) 

The analogous vector bp^((3] can be calculated in the reference frame of |3. The 
two vectors are related by the formula 

b^p(a)=U^(a,(3)b^J(3) (4.16) 

By construction these vectors are orthogonal to the face f^p so that, from a math- 
ematical point of view, they live on the dual Voronoi edge (a, |3) (see appendix 
0). We choose the order of the the edges vectors E^^ in formula |4.15| in such a 
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way that b^p((x) points toward the outside of the simplex oc. Among the vectors 
b^pla) of the same n-dimensional simplex there is the following identity 

n+l 

^b^p(a)=0 . (4.17) 

The antisymmetric tensor V'^' a, ailoc) can be written as a bivector of the 

b^p((x) 



(4.18) 

where oc — 1 and a. + 1 indicate the two n-simplices which have a face and the 
hinge h, in common with a. The volume V(a) of the simplex a can be written as 
function of the b^p in the following way 

(V((x))--^ =^e,,..,,e'--''^Sb,^U(3)...b°:^((3) , (4.19) 
where the indices (ji, jn) run over all possible values of the dual Voronoi edges 



(a, (3) ,(3 = 1,..., n+l and j is a fixed index. The equation |4.19| does not depend 
from the choice of the index j in the e-indices, or, in other words, it does not 
depend from the b that is left out. 

As said in the previous section we can put the origin of the reference frame of 
a in the circumcenter . If ( a) , a = 1 , . . . , n, i = 1 , . . . , n + 1 are the circumcentric 



coordinates |]8§] of the n + 1 vertices of oc we have 

n+l 

^z^ioi]=0 . (4.20) 

In these coordinates we can express, as it is easy to see, the vector b?(a), i = 
1 , n in the following way 

This equation can be inverted, so we have 



1 



(n+1)!(n!V(a])^ 



■a, . 



--On- 
t-ik 



(X 



.bf-- 

Jn — 



oc 



(4.22) 
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These two equations |4.21| |4.22| show that there is a one to one correspondence 
between the circumcentric coordinates of the vertices of a simplex a, and the 
bi(a). 



The gravitational action |0| can be written, by the equation |4.18| , as a function 
of the bf (a). Anyway we are still dealing with a second order formalism, since the 
connection matrices U^[(x,, (3) and the bocp((x) are both functions of the coordinates 
of the edge of the simplicial complex. Now we propose a first order formalism in 
which U{J(a, (3) and botp(a) are independent variables defined on the Voronoi edge 
(a, (3]. We consider equation |4.16| as a constraint of the theory. This constraint 
fix only n degree of freedom for each Voronoi edge (a, (3), and not the ' 
conditions that are necessary to determine Ug((x|3). Pairwise |4.17| is the second 
constraint of the theory for each Voronoi vertex. 

In the first order formalism we consider connection matrices more general than 



the Levi-Civita matrices defined by equation |0 



This implies that in these cases the gravitational action in the form ( |4.8|) could 
be dependent from the reference frame in which it is written since it is no longer 
true an equality like |4.13| . So if we define the following antisymmetric tensor on 
the plaquette H as 
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+ U,p...U5H 5HVf^'(6l:)U6H6H Up„) . (4.23) 

K — Ik kic — I 'Cl C2 

the action can be written in the form 

S ^ --^Tr [u'^^W'^ioc]) (4.24) 

The above action coincide with the previous action in second order formalism 

and 

since (in matrix notation) 

wC^Hcx) = U,,p...U5,_,5,W'^'(50U6,5,_,...Up^ (4.25) 

the action ( |4.24|) result to be independent form the frame in which it is written 
essentially for the same reasonings of the lemma (|4.11|) . 

Moreover the action (|4.24|) is invariant under the following set of transformations 

Uccp ^ u:,p = 0(cx)U,pO-^(|3) 
bap(a) ^ b'^^ioc) = 0{oi)\}ocfi[oc] (4.26) 

where 0(a) and 0(|3) are two element of SO(n). 

We want to close this section by adding some consideration on the nature of 
the b^p. In the interior of every simplex ct, as we have remarked above, the space 
is flat. So in every point of the cotangent space to a the metric is flat 

g(ff) =6abe"(o-)®e^o-) (4.27) 

where e'^(a),a = 1,...,n form an orthonormal dual base in the cotangent space 
of ff. Of course it is evident that in every point of the cotangent space of cr we 
can choose the same base e''(ff). In the form (|4.27|) the metric coeflicients are 
automatically diagonal and so we can say that trivially a n-dimensional simplex 
is a local inertial reference frame. 

It seems more natural to choose as base the edges of the simplex cr because this 
makes contact with the geometry of the simplex cr. Let 



E^[a) [1= 1,...,n 



(4.28) 
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n linear independent edge vectors. The other edge vectors can be expressed as 
linear combination of the n vectors above. Associated to these n-vectors we can 
fix (uniquely) n covectors of the dual base such that, by the definition (|4.15| ), 

b^(cT)(E,(cT))=5i^ . (4.29) 

where 

b^((T] = (4.30) 
rL!V(0j 

The transformation law between the couple of the orthonormal base ea[<y) and 
their respective duals e''(o-) and the couple of E^(o-) and h^[(r) are given once we 
fix the component of ET,(cr) and h^{a) in the base ea(cr) and e'^(cr), more precisely 

e"(0)=E^(0)b^(o-) 
eja) =b^((T)E^((T) . (4.31) 

In these basis we can write the metric tensor as 



g(o-)=6-Sti(o-)b^(a)E^(cT)®E,(cT) 

g(a) = 5abE^(a)E^(a)b^^(a) ® h^{a) (4.32) 

It is now clear that we can do the following identification respect to the continuum 
theory 



e^(x)^b^(cT) , (4.33) 

where with e^(x) and e^(x) we indicate the n-bein of the continuum theory. The 

discrete n-bein, in analogy to the continuum case, are determined modulo the 

action of the orthogonal group because if we perform an orthogonal transformation 

on E^(a) and b[;J^(cr), the equations ( |4.32| ) remain form invariant. 

Moreover the second equation of ( |4.32|) says that the coefficient of the metric 

tensor in the base of the b*^(cr) are the usual metric coefficients written in Regge 

calculus 



(4.34) 
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where, in the standard notation, = |E^P and with = E^ — E^, l^^ = |E|j;vF- 

4.4 Remark on the Orientation 



In the definitions ( |4.18|) , ( |4.23|) , (|4.5|) we have tacitly assumed that the boundary 



of the plaquettes is oriented. As a consequence, for example, we can establish that 
the Voronoi-edge (a, (3), in the direction from a to |3, is positive oriented, while 
the opposite Voronoi-edge |3(x is negative oriented. So this naturally fix a rule to 
define the bivectors V*^' '^^'^^[a], W^[a) and the holonomy matrix Ujj^ along the 
positive direction. In particular we use the notation V'^' '^^'^^((x) = V''^' ap^l*^) 
stress that in the definition |4.18| the positive direction is from a to (3 in order that 
in the bivector we have to use first the vector bccplo.) and after bcc6|^(a). Pairwise 
we indicate as W+ the antisymmetric tensor |4.23| that in these new notations is 



+ U«p...U,K_^,HVf^;^(6]:)U,H,H_^...Up«) , (4.35) 



and U+^^ = U.^^, that is the holonomy matrix in ^ is defined along the positive 
direction (a, (3). 

The negative direction is from a. to 6]^, and as in the definitions above we have 



and 



+ U^H...U5HpV^^^(P]Up5H...U5Hj^^^^ , (4.36) 



l-l aa = ^a5]2l-^6]25h_T •••U(3a • (4.37) 

Now we are going to prove the following lemma: 

Lemma 4.4.1 The action \4-24 is independent from the chosen orientation 
for the boundary of the plaquettes H in the sense that 

S^-l^Tr (u+lw+^cx)) =-l^Tr (u-^W-^a)) . (4.38) 

K K 
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Dim: since the trace of a matrix and its transposed are the same we can write 



(4.39) 

Now we have that 

U+L"' = U,5hU5H5H^...Upo, (4.40) 

and 



+ Ua|3...U5h_^6hV*H6h_^(6]^)U5h5h_^...U|3, 

Consider the orthogonal matrix 



r — Ua(3U|35h...U5h §h 



and manipulate the last trace in 4.39 in the following way 



(4.41) 



(4.42) 



1 

2 



A straightforward calculation shows that 



(4.43) 



1 



+ U6H5h_^...U|3aVS(0^)UaP-U5H_^«H 



and 



.u 



a5^ 



(4.44) 



(4.45) 



^--i — 0.0. ■■ I '2_~~ \~' °k°k ' ' ^-k-)) (4.46) 

Since the trace |4.46| is independent from the reference frame where it is written, 
we have the equality |4.38| . 



These last two equations imply that 
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This property of the action is necessary not to have ambiguity in defining it. 
As we have seen on each hinge there is no preferred orientation (c.f section ^) 
and, as is easy to see, no orientation on the corresponding Voronoi'-dual plaquettes 
even if the manifold is oriented. 



4.5 First Order Field Equations for Small Deficit 
Angles 

As remarked in [^, in the second order formalism the action (|4^) is equivalent 



to the Regge action for small deficit angles 9h (since in this case sinQh ~ 9h). In 
the first order formalism we don't have angles 9h, but the only variables related to 
the deficit angles are the connection matrices U^^p. Then we assume, by definition, 
that the small deficit angles approximation in the first order formalism is the 
passage from the group variables Uap of SO(n) to the algebra variables (^ae> of 
so(n). As a consequence the connection matrices can be written in the form 



Uccp = I + e ct)ccp + o(e) (4.47) 

This implies, as it is easy to verify, that the only gauge transformations ( |4.26|) 
which can be compatible with the equation (|4.47| ) are of the type 0(a) = I, or, in 
other words, the approximation in the form ( |4.47|) is also a gauge fixing. 
In order to avoid technical complication, that we shall discuss in the next section, 
we now suppose to substitute the constraint bccplcx) — U(xpb|3ct((3) in the action 
for each Voronoi edge ( cx, (3 ) . 

A straightforward expansion of the action up to the first order shows that 

S = -le }^Tr((cl)„p + cl)p5H... + ct)6Hj°W^a))+o(e] (4.48) 

where '^^ is the bivector '^^ tj^g gero order in which we have done 

for each matrix the approximation (|4. 4 71) and for each Voronoi'-edge we have solved 



the constraint (|4.15| ) up to the first order, that is to say 

bap(a) = (I + e ct)<^p)bpj(3) + o(e) (4.49) 

Demanding that the action be stationary under the variation respect to (cj^ixp) we 
obtain 
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= e y_ °W^^i^^ +o(e) =0 (4.50) 

where he (a|3) means the sum over all plaquettes h to which the Voronoi'-edge 
a(3 belongs. The tensor is an antisymmetric tensor of type (2, 0] in 

n-dimensions, so has ^^^^^^ independent components, the same independent com- 
ponents of 4)^13 that have to be determined. 

Let's consider the face between the simplices a, and (3, dual to the Voronoi'-edge 
a|3. The equations ([|72[) that determine the Levi-Civita connection between 
the reference frames of the two simplices , to the first order (14.471 ) 



{z'^_,[oc)-z^[c^)) = (l + ecl),p)^(zti(P)-z|;(P)) . (4.51) 

The zf ( a) , a = 1 , . . . , n are the coordinates of the vertices i = 1 , . . . , n of the face 
in the reference frame of the circumcenter in a. The circumcentric coordinates 



4.20| of n + 1 vertices of a n-dimensional simplex are in one to one correspondence 
with the n + 1 bf((x), {i = l,...,rL+ 1} . Equations ( |4.51| ) can be seen as an 
equation for determining 4)ccp as function of the zf (a) (and then of the (a)). 
Each cjjccp, being an antisymmetric matrix in n-dimension, has "^^^^^^ degree of 
freedom. The number of independent equations are n(n — 1). Anyway there is 
the following identity among the edge components of the n — 1 -dimensional face 
(a,|3) 

Y_ V^^"'Hoi) = , (4.52) 



The (14:521) fixes ^^^^ degrees of freedom, so that the (^^ ) together has 



^ 2 linear independent equations. By linearity there is for cj^ccp only one solution. 
Anyway this is the connection that we adopt in the second order formalism. We 
want to see if this connection, that is the Levi-Civita or Regge connection, satisfies 
equation (|4.50|) . If the connection is Regge we have (in matrix notation) 



V'^'(a] = U„pV'^ni3)Up,= ... = U,p...U5K^5HV"^H5]:)U5K5H_^...Up, (4.53) 
so that to the zero order 
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V'^V) =V'^'(|3) + 0(e) = ... = V(^H5]:) + 0(e) . (4.54) 
These facts imply that 



e Y_ °W^^"''= Y_ (eV'""'^Mcc) + eO(e)) =0 + O(e^) (4.55) 

he(aP) He(a(3) 

that is to say that Regge connection is solution of our first order equations for 
the connection matrices in the limit of small deficit angles. 



4.6 First Order Field Equations: the General Case 

In the previous section we have seen that in the case of small deficit angles the 
Regge calculus is solution of the first order field equations. 

We are now going to deal with the general problem. We want to derive the 
equation of motion by varying the action respect to Uccp and b^p. This formula- 
tion of gravity on lattice is close to the popular Ashtekar's variables [|6| in non 



perturbative canonical continuum quantum gravity (see also [|3 for a covariant 
version). Of course we have to take in account the constraints (|4.16| ) and (|4.17|) , 
so in order to perform independent variation of U^p and b^p it is necessary to 
introduce the constraints in the action by using Lagrange multipliers. Then the 
first order action will be 



S = -^}^Tr(u^^W^cx))+}^A,p(b,p(a]-U,pbp„(l3]) 

+ }^Tr(A(«P'(U„pUTp-I))+^^(a] (^b«p(a)) (4.56) 

(ap] a \p=i / 

where A''^'^' and |j.(a) are n-dimensional vectors and are Lagrange multipliers. 
^(aP) is an n X n matrix. The constraint UccpU^^p — I is introduced to restrict the 
variation of Uocp on the group SO(n). 
Let's introduce the following lemma: 
Lemma: The action, in matrix notation, 



S' = Tr(AA)+Tr(A(AAT-I)) (4.57) 
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gives the following equation of motion if we assume that the variation 
respect to A be stationary 

(AA) = (AA)T . (4.58) 

Dim: If we consider the variation of the action respect to A and A^, using the 
property that Tr(M) = Tr(M^), we have 



6S' = ^Tr(5AA + A"^6A"^) + ^TrA(6AA^ + A5A"^) + ^Tr(6AA"^ + A6A"^)A"^ (4.59) 
that implies 

= aa + a + at = o 

= (AA)"^ + A^ + A - , (4.60) 



(5A)AT 
A(6AT) 



from which, subtracting term by term these two equations we have (|4.58|) . 

We can apply the lemma to the action (|4.56| ) for the variation respect to Uocp, we 

obtain the following field equations 



Y_ (U^«W^a))tj-A,p,b,p(a)^= Y_ (U^ocW^(<x))X - A,p jb,p((x) , (4.61) 

H6(a(3) He(a(3) 

The next step will consist of determining the field equations for the variations 
of b^p. For this aim it is necessary to determine the quantity We recall 

that in the circumcentric coordinates, the coordinates of the i-th vertex ((3) is 
given by formula |4.68| 

We can write the formula |4.19| in a way not dependent from the chosen index 
j. If we call (V'''(cx))^^^ equation [4.19| , we have that it can be rewritten as 



(4.62) 



in such a way that it does not depend from any index. Evaluating the derivative 
of this last expression we get. 



V(<x)) 



1 



3bf ( a) 



n!(n-l)V2(cc) 



(4.63) 
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We are now ready to derive the field equations for the variations of the bap. 
The starting action is (|4.56|) . Furthermore we have remarked that the action is 



invariant under the orientation of the boundary of each plaquette, so we can write 
it in an orientation independent way 



H (a(3) 

+ _^^(^x) (^b,p(a)] . (4.64) 



v|3=1 

A straightforward calculation show that: 
with the notation 



Ka - = (4.65) 



+ At ap-^,((x) =0 . (4.66) 



In the field equations (|4.66|) and (|4.61|) we have to determine also the Lagrange 



multipliers A^p and \x)-[oi). We now propose a method for determining them. In 
each simplex a, or, equivalently, in each point of the dual complex, we have n + 1 



vectors b^p(a),b^i (a), ...,b^n((x) and among them there is the constraint (|4.17|) . 
Let us choose n of this vectors for example b^|3(a), ...,b]^5^ Ja], in other words 
we are solving the constraint (' |4.17|) . This is a base in the tangent space of the n- 
simplex cx so that we can write each Lagrange multiplier, such that Ajj^p, considered 
as a vector in tangent space of a as 

A^p = c>^p(a)...c-pbX_^ (a) (4.67) 

The set of n vectors b^p((x), ...,b^^^ (a.) are in correspondence with the cir- 
cumcentric coordinate 



1 



(n+l]!(n!V(a)]^ 



-a, ...an_i 



- On-l T^ai 
k "^i, 



(4.68) 



4.6 First Order Field Equations: the General Case 



55 



where the indices (ji, jn-i,T<) run over (x|3, a6i, aS^-i and in this case n+l = 

From a geometrical point of view elementary arguments show that z^_^-^[^x) is the 
circumcentric coordinate of the vertex of a which is the intersection of the n 
n— 1 -dimensional faces whose normals are b)^p(a), ...,b)^5^^(a). 
Pairwise to each bocsjoc) we can associate a circumcentric coordinate z^. (a) ob- 
tained as in the definition |4.68| in which n + 1 i-^' i = aSt and (ji, jn-i, !<-) now 
run over a|3, (x6i, (x&n in which, of course, ocSi is omitted. In this way we can 
make a one to one correspondence within each 'bo^&^[o^) and the edge of the simplex 
a Eccsjoc) defined in the circumcentric coordinate 



ES5j<x)=z^Ja)-z:^+,(a) (4.69) 
It is straightforward to see that 

(b„6j(x))"(E^,(cx))^ = n!V(a)6y . (4.70) 

Consider the equation (|4.671 ), project it along the edge Ejj^p and repeat the 
procedure for each Ejj^g. ,i=l,...,n— 1. We will obtain 



1 

n!V(a) 
1 



In other words for determining Accp we have to determine its projections along 
the above edges. For this purpose we project equation ( |4.61| ) along E^g. (a) and 
E'^l3((x) and we find that 



A^pE„5, H(a) = --^ Y_ ((U^„W^a))pq-(U^,W^(x)),p)E^Ja)E^p(a) , 

(4.72) 

which fixes the coefficients c^p, i = 2,...,n. A problem arises with the coeffi- 
cient cj^p, in fact equation (|4.61| ) does not determine the component of A^p in the 
direction of bj^p(a), because the following transformation A^p i— > A^p + (xbj^p does 
not affect the equation. It seems necessary to use the field equation for b^p for 
evaluating c^p. 

Anyway if we project (|4.61] ) only along E'^p(cx) we obtain that 
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(4.73) 

which, substituted in the starting equation, gives the field equation for Uocp, 
without Lagrange multipliers 

xEl^[oc]K^[oc) - Y_ ((ULW^a))3k-(U^,W^a)]^j)E^p(a)b^,p((x) (4.74) 

He(a|3) 

Of course the reasoning for the Voronoi-edge a(3 can be applied to other Voronoi 
edges (xbi. The unknown coefiicients will be c^j , 1 = 2, ...,n— 1. 
In the reference frame ex. we write the equation for the variation respect to bocp 
and bocSi in the following manner 

Ks.+Ks. + ^'^ioc] = (4.75) 

where in the symbols B we have summarized all the terms in the equations ( |4.64|) 
which don't contain the Lagrange multipliers. We recall that the expansion of A^p 
in the base of b is given by (|4.67|) , pairwise 

= cU^ap + - + cS^as, + - + , (4.76) 

Subtracting the two equations ( |4.75|) between them, we have 

(B^p-B^5j + (c^p-c^)bSp + ... + (c-p-c«5„_,)bS5,_, ^0 . (4.77) 
Projecting along (Ec^plajja we get 

^'^^ = ^'^'^ ~ ^^r^^^^P " B^5j(E,p(a))a (4.78) 
in this way we determine c^p, and a similar procedure gives 

''^^^ ^ ^"^^ ~ ^^iVR^^^^ " B^p)(E«5ja))a (4.79) 
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D 




C 



A B 

Figure 4.2: The Most Elementary Triangulation of the 2-Sfere 

Once we have determined A^p and A^. we can determine \i.°^{(x.] by summing 
the two equations 

^t'^(cx) = ~{K^ + ASsJ - (B^p + BS,J (4.80) 

We have used the equations for b^p and bc^^ to determine A^p, ^abi, and 
To determine b^p and ba&i we can solve the equations for bpa((3) and b5^oc(5i) in 
which the Lagrange multipliers |J.'^(|3), ia.'^(6i) have been determined by two couple 
of different equations. Then we will use the constraint baploc) = UccpbpccO) and 
'^ocSiM = ^aSi^Sioci^i) ^^d SO on. In this method, at the end, the only unknown 
quantity is the Lagrange multiplier A^^ corresponding to the Voronoi'-edge (x6n 
and to b<x6n- Anyway we have that bccs^ is determined once we know all the other 
b and we can use the equation for b^t^^ to determine \a&n function of the other 
b. 

Let do an example that show how this method work in practice. Consider the 
boundary of a tetrahedron as a two dimensional triangulation of the sphere 

Let call a the triangle ABC, ^ the triangle BCD, y the triangle ACD and 6 
the triangle ABD. 

The edge BD is the face (one-dimensional) between a and |3, bapla) will be the 
vector in a perpendicular to this face. DC is the face between (3 and y, and bpydS] 
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the vector perpendicular to it. AD is the face between y and a, and bycc(y) is the 
vector perpendicular to it, and so on. AC is perpendicular to the Voronoi-edge 
ay, BC to (|35), AB to ocS. 

Let's write the equation of motion for the variation of the b for this system: 







= 


edge 


((x|3) 
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edge 


(ay) 
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edge 
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Bpa- 
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edge 
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edge 
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edge 




ByP - 
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(y(3) 




(A«^U,^]"+^'^(y] 
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edge 


(ya) 




+ + ^t'^(y) 


= 


edge 


(y5) 






= 


edge 


(5a) 




-(ApsUpsl'^+^^lS) 


= 


edge 


(5P) 






= 


edge 


(5y) 



Prom the first two equations we choose to determine A^p, A^^ and (J.'^(a) by the 
above method. The equation for the Voronoi-edge a6 will be used to determine 
A^5. The equations for the Voronoi-edges |3y and (36 will be used to determine Ap^, 
Aps and |J.'^((3). The equation for the Voronoi-edge |3a will be used to determine 
bpcc- Prom the constraint b^p — Uapbpcc we determine bocp too. The equation for 
the Voronoi'-edge y(3 may be used to determine (J.'^(y). The equation for ya may 
be seen as an equation for determining b^cc, and so we determine bccy too. The 
equation for y6 may be used to determine Ays. The equation for 5a determine 
|j.'^(5), those for 5|3 and 5y determine bsp and bsy and so also bps and b^s. The 
other b are determined from the constraint on the b ( |4.17|) . 

In this way we determine all the Lagrange multiplier as functions of the Uap 
and bccp. Using also the constraint we have that all the field equations fix a number 
of degrees of freedom that is equal to the number of degrees of freedom of Uc^p 
and bap . 
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In this section we shall discuss the quantum measure to associate to the previous 



classical action. The action is invariant under the action of the group SO(n), 

the gauge group (|4.26| ). 

Let us fix the following notation: we define 



(4.82) 



and with 



^(Uc,p) (4.83) 
the Haar measure of SO(n). The partition function for this theory will be: 



n+l 



(4.84) 



(3=1 



f^S (bap(a) -U„pbpa((3)) M.(U„p)M,(bcc(3] 

a(3 



in which the first product fla a product over all vertices of the dual complex, 
and the second product Y[ae> product over all the edges starting from the 

vertex oc. The second delta function 6 (bcc(3(a) — Ucc(3b|3oc((3)] has to be introduced 
one time for each edge, so that if we have introduced it for the edge (x|3, we have 
not to introduce it for the edge (3 a in the reference frame (3. It is straightforward 
to see that the measure of the partition function is invariant under the gauge 
transformation ( |4.26|) . In fact if we perform a gauge transformation of the type 
( |4.26|) , the modified terms of the measure are: 



'n+l 



(4.85) 



^P=i 

Now |J.(b^p) is equal to det (0(cx)] ^(bccp] and so is equal to \<i{hot^]. The Haar 
measure of SO(rL] is right and left invariant so M-lU^p) — \x{U.otp,]. By the equations 
( |4.26|) , we have that the two delta of (|4.85|) can be written in the following form 



TL+l 



5 0(a)(}^b^p(cx) 



(3=1 



6 0(a) bocp(a)-U«pbp„((3) 



(4.86) 
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that, from the properties of the delta function, will reduce to 1 over the determi- 
nant of 0(a) for the deltas of the quantities in parenthesis that multiply 0((x). 
Then we have that equation (|4.85|) can be written as 



S \ X. ^ocfi[oc]j 6 (bap((x) - Uo,pbp J(3)) ^(U<^p)^(bocp) , (4.87) 
and then the invariance of the measure under gauge transformations. 

4.8 Coupling with Matter 

In the continuum theory on Riemannian manifolds the coupling with fermionic 
matter is given by the following term of the Lagrangian density 

£ = (^\>e^'Y^V^^\> - V^e^iJ^y>) + mi|;ip (4.88) 
where y*^, a = 1 , ...n are the Dirac-matrices satisfying the Clifford algebra 

yV+y^y" = 26''^ (4.89) 

i|) the n-dimensional Dirac spinor field (if) = the covariant derivative, 

the n-beins that is, on the tangent space of the Riemannian manifolds (M, g) 
where the Lagrangian density (|4.88|) is defined, the vector fields such that 



g^^-(x) - ett(x)e^(x)5-^ . (4.90) 

We are assuming that the Riemannian manifolds (M, g] in question have a spin 
structure, that is the second Stiefel- Whitney class is zero. 

We have now all the ingredients to define the coupling of gravity with fermionic 
matter on the lattice in analogy with continuum case (for related works see 
H S). 



Let "V = 2rL or T = In + 1 (depending if n is even or odd) and consider the 
2^ dimensional two-fold covering group of SO(n]. So instead to consider the 
connection matrices Uap, we have to consider the 2^ x 2^ connection matrices 
Dcc|3 such that 

Dapy'^D^J = (U,,p]^y^ (4.91) 
where in the left hand part of ([4.91j ) we have used matrix notation. Given Dc^p we 



determine 11^^, but if we know we determine D^p up to a sign. In particular, 
as it is easy to derive from (|4.91|) , we can write Uc^p as (0) 
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(Uap)ab = ^Tr(yaDapyb) (4.92) 

In the discrete theory we assume that the spinor field is a 2^ complex vector 
defined in each vertices of the dual simplicial complex, that is to say a map that 
to each vertex oc associates the 2^ complex vector iKa). 

In order to define the covariant derivative on lattice we have to derive the distance 
|(X|3| between the two neighborly circumcenters in a and |3. The distance AHi of 
the circumcenter in a from the face a|3 can be calculated by deriving the volume 
of the TL-dimensional simplex obtained by joining the circumcenter with the it- 
vertices of the face a(3 and dividing it by n and the volume of the face itself. So 
we have that 

AH, = -l ^^',^°^it°'"' (4.93) 

in which zf are the circumcentric coordinates of the vertices of the face a|3 and 
with |bap(a)| the module of the vector that divided by (n — 1)! is equal to the 
volume of the face. 
In the same manner we have 

AH. = J,&^M^ (4.94) 

At the end we have |a|3| = AHi + Ahz- 

At this point we are ready to define the covariant derivative (V^ip)(a) on 
lattice 

(V,^)(a)^ °'^^'PJp-1""' . (4.96) 

So far the discrete version of the the action for the coupling between gravity 
and fermionic matter can be written in the form 



oc (aP),p=1,...,n+1 ''^l ' 

- D^^M^)hle.yaMoc])+mMoc]Moc)) (4.96) 
so that the quantum measure that includes also fermionic matter can be written 

as 
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(S+Sf 



a|3 



(4.97) 



where S is the action for pure gravity ( |4.24| ), and |j.(Dccp) the Haar measure on the 
two- fold covering group of SO(n) while (j.(i|i(cx)) = dip (a) the standard measure 
on . 

The coupling with Fermionic matter on lattice presented in this paragraph is 
quite similar to the one introduced in reference [|8^. The differences consist of the 
fact that the action is here group theoretical and not Regge, and this formalism 
implies that from the beginnings the coupling is written on the dual Voronoi 
complex. 



Appendix A 
Branched Polymers 



The notion of Branched Polymers is crucial in the classification of the so called 
" Elongated Phase" of 4-D simplicial quantum gravity. We are now going to review 



5§] some basic concepts of statistical mechanics on graphs. 

An abstract graph is a set of points called "vertices" connected by lines. These 
lines introduce a natural concept of connectedness. A graph can be seen as a one 
complex. Prom now on all our graphs will be abstract, not embedded in some 
affine space so that we may omit the adjective abstract. A graph is simple if two 
distinct vertices are connected at most by only one line and the vertices at the end 
of each line are always distinct. A graph that is connected and simple is called a 
tree. The number of lines 1(1] at the vertex i is called "the coordination number" 
of the vertex i. In a tree with n vertices we have that Y.i=i Ui-) = 2(n — 1 ). A tree 
is rooted if we label one vertex which is called the root. 

Two trees T and T' are combinatorial equivalent if there is a one to one map 4> 
between them such that the vertices of T are mapped uniquely into the vertices 
of T' so that a line joining two vertices of T', v( and v|, is a line of T' if and only 
if it is the image of a line of T by ({j, in the sense that there exist two vertices of 
T V| and vj joined by a line of T whose images by cj) are v[ and v- and the line 
joining them is mapped into the line joining v( and vj. If the two trees are rooted 
they are combinatorially equivalent if there is a one to one map such that together 
with the properties described before it maps the root of T into the root of T' and 
viceversa. In statistical mechanics tree graphs are called Branched Polymers. 
The statistical mechanics of these abstract tree graphs is also known as mean 
field theory of Branched Polymers to distinguish it from the case of embedded 
Branched Polymers. 

Let's define as r*(n) the number of combinatorial inequivalent rooted tree 
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graphs with n vertices in which the root has coordination number equal to one 
and all others vertices have coordination number less than or equal a. By Caley's 
formula we have 



r«fn) 



(n-1) J-Jdt-l)!™ (n-1)! " ^ 



i=2 

Let £,'^(n) indicates the number of combinatorial inequivalent tree graphs with 
n vertices whose coordination number is less than or equal to a. We have that 

Prom the above formula or from geometrical arguments it is easy to see that 

l°'{yv) = -r°'[Jy + ^] . (A.3) 
n 

The partition function R*(|3) over the ensemble of rooted tree graphs, whose 
vertices have coordination number less than or equal a, is defined in the following 
way 

oo 

R-^dS) = |3^-V«(n) , (A.4) 

n=2 

and the analogous 2*^(13) unrooted 

00 

Z'^((3) = _^(3^r'^(n) . (A.5) 

TL=1 

Prom these definitions and the identity [A.3| we have that 

^z-(3) = ^R-(3) . (A.6) 

As regard the definition of the Green function for tree graphs with coordination 
numbers up to a, the strategy is completely similar to the same case for dynamical 
triangulations (c.f. section |3.51|) . A path between two vertices of a tree graph is 
given by a 

sequence {\]X=\ of vertices such that the vertices j and j + 1 (j + 1 7^ j and 
j — 1,...,l— 1) are joined by a line and ii coincides with the first vertex and ii 
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with the second vertex. The length of the path is the number of its line, I — 1 
in this specific example. The distance between two vertices of a tree graph is the 
length of the minimal path that joins them. 

As usual the micro-canonical Green function is the sum over all combinatorial 
inequivalent tree graphs T'^(r, n] whose vertices have coordination numbers up to 
(X with n vertices and with two labelled vertex at distance r 

G'^lr.n)^ Y_ (A.7) 

Ta(r,n) 

and the relative grand-canonical partition function is 

oo 

G-(r,(3)^}^rG°^(r,n) . (A.8) 

TL=1 

The relative susceptibility function is defined as follows 

oo 

X'^((3] = }^G«(r,|3) (A.9) 

T=0 

As for the dynamical triangulations in the definition above we can exchange 
the sum over r with the sum over n so that 

T T°'{r,rL) 

is the sum of all combinatorial inequivalent tree graphs with coordination num- 
ber up to a and with two labelled vertices. Since for large n asymmetric tree 
graphs prevails this sum is asymptotical equal to n^G'^(r,rL), so that as in p.29| 
the susceptibility function is asymptotically equal to 

X"(13] X ^Z-(13] . (A.ll) 
Now we are going to use the following recursion relation 

Lemma A. 0.1 The function r"(rL) can be written as 

oc-l , y 



y. 

Y=0 1^1 ,...,ny>2 1=1 
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Branched Polymers 



This recursion formula can be proved straightforwardly or by using formula |A.l 
or by geometric considerations. Anyway from this lemma we have the following 
evident corollary 



Corollary A. 0.1 



R-((3) = |3 V -[R-(P)] 

Z -v/l 



y=0 



y! 



(A.13) 



At this point [A.13| can be used to derive the differential equation for R'^IP), 



(A.14) 



This equation is divergent at the point |3'^ where the denominator vanish, that is 
to say 



cc-1 



(A.15) 



This equation has only one solution for each value of ex. Since is divergent 

in (3'^ we have that 



dR'^ 







(A.16) 



By the implicit function theorem (3(R'^) result to be analytic around R"(|3'^) so 
that we can expand it in Taylor series 



(3(R'^) - (3^ = ^^^^(R'"- R''(|3'))^ + o ((R'^- R'^(|3^))^ 
and for (3 near |3'^ we get 



(A. 17) 



R'^lpjxR^lp^j + Clp^-p)^ (A.18) 

These considerations together with equations |A.6| and [A.ll| imply that for |3 
near (3'^ 



X"((3)x((3^-P)-I 



(A.19) 



so that the susceptibility of Branched Polymers has the same exponent ys = ^ of 
the mean field. 



Appendix B 
Baby Universes 



The concept of Baby Universe in dynamical triangulations has been introduced for 
the first time by S. Jain and S.D.Mathur The original idea of these authors 
has been to introduce a way for describing the roughness of triangulated surfaces. 
Successively the average number of Baby Universes has been used as a sort of or- 
der parameter to classify the two phases, crumpled and branched polymer phases 
which have been seen for the first time in the Monte Carlo simulations of 4-D dy- 
namical triangulations. In fact in the articles and |j5^ numerical simulations 
have shown very few baby universes in the crumpled phase and a cascade of baby 
universes in the branched polymer phase. Motivated by these data it has been 
tried to put in correspondence the baby universes with the branched polymers. 
The deep reason of this correlation is discussed in chapter ||. Anyway seminal ar- 
guments in favor of the connection between branched polymers and baby universes 
of minimal size can be found in the article [|8| . 

Let's give the main definitions regarding Baby Universes in the spirit of the 



article We shall restrict, for convenience, to the four dimensions and to the 
case in which the baby universe has the neck of minimum size (Minbu). A Minbu 
can be considered as a triangulation of a PL manifold in which we can distinguish 
two parts: the Mother universe that is the part that has the majority of the volume 
of the triangulation, and the baby universe which is the part with the minority of 
the volume. The two parts are glued together by the boundary of a four-simplex 
that is called the neck of the baby universe. 



Now we are going to introduce, following a way to enumerate the average 
number of baby universes with given volume V in the canonical ensemble of all 
triangulations of a PL-manifold with N4 four-dimensional simplices. Anyway we 
want to stress that this enumeration is not rigorous but it is in good agreement 
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Baby Universes 



with the numerical results [|5^, [^. Suppose that the topology is that one of 
the sphere S"^ since our future discussions will be concentrated on this case. Imag- 
ine to separate the mother universe and the baby. We will obtain two triangula- 
tions of the sphere S"^ with a boundary made by the boundary of a four-simplex. 
The mother has TM4 — V four dimensional simplices and N2 two-dimensional sim- 
plices, the baby V four-dimensional simplices and N2 two-dimensional simplices. 
Obviously the sum N2 + N2 will be equal to N2 + 10, where N2 is the number 
of two-dimensional simplices of the original triangulation, because in separating 
the mother and the baby we are counting the bone of the neck twice. Suppose 
now to fill the hole of the boundaries by gluing a four simplex to each of them. 
So that the number of combinatorial inequivalent triangulations of the mother 
universe is now PalS"^, N4 — V + 1,N2) (c.f. |2T^ ) and that one of the baby is 
Pq(S'^, V+ 1 , N2). The antsaz is that all the combinatorial inequivalent configura- 
tion with a Minbu are obtained by marking a four simplex in the mother universe 
and one in the baby, and by gluing them through identifying the two marked 
simplices. We suppose, as usual, that we are in the regime of large N4 so that 
asymmetric triangulations prevail. This implies that the marking process will pro- 
duce (N4— V+1 JpqIS"^, N4— V+1 , N2) combinatorial inequivalent triangulations for 
the mother universe and (V+ ])pa{S'^,V + 1 , N2) for the baby universes. So that, 
at the end, the number of the combinatorial inequivalent Minbu M(S'^, N4, N2, V) 
is 



M(S4,N4,N2,V] = 5!(N4-V+1)(V+1) (B.l) 

X pJS^ N4 - V + 1 , N^)pa(S^ V + 1 , N^') . 

where 5! comes from the fact that in the case of asymmetrical triangulations there 
are 5! combinatorial inequivalent ways of joining two four-simplices. Any Minbu 
of the characteristics defined above can be obtained in this way, so that |B.2| really 
can be considered as a formula that gives an estimate of the number of all Minbus 
of baby-size V and number of bone N2 in the micro-canonical ensemble of all tri- 
angulations of S'^ with N4 four-simplices and N2 bones. Consider a triangulation 
T of S"* with N4 four-simplices and N2 bones, which has a Minbu of size V and N2 
bones. Let us indicate by nTlS"*, N4, N2, V) the number of combinatorial inequiv- 
alent triangulations that can be obtained from a given triangulation T by cutting 
a baby universe and attach it on another four simplex of T, previously marked, as 
described above. It is clear that M(S4, N4, N2, V) = Lj'n-TtS'^, N4, N2, V), so that 
the average number of combinatorial inequivalent Minbu in the micro-canonical en- 
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samble can be approximated by dividing M.(S'*, N4, N2, V) by the micro-canonical 
partition function PalS"^, N4, N2) 



Ar(S4,N4,N2,V) ^ 5!(N4- V+1)(V+1) (B.2) 

Pa(S^ N4 - V + 1 , N^jpJS^ V + 1 , N'^] 

Pa(S4,N4,N2) 

In the canonical ensemble the average number of Minbu with volume V can be 
calculated in the following way 



Ar(S4,N4,k2,V) ^5!(N4-V+1)(V+1) (B.3) 
Ln, e''^^^ Pa(S"> N4 - V + 1 , N^)p JS^ V + 1 , ] 

Z(S4,N2,k2] 

which using the relation N2 + N2=N2+10is equal to 



A/'(S^N4,k2,V] f^5!(N4-V+l)(V+1)e-i°'^2 (B_4) 
Ln^ e*^^"^^ Pa(S^ N4 - V + 1 , N^) e^^^2p^[S\ V + 1 , ) 
' Z(S4,N2,k2)' ■ 

Finally 



A^(S^N4,k2,V) ^ 5!(N4-V+1)(V+l)e-^°'^^ (B.5) 

Z(S4, N4 - V + 1 , k2]Z(S4, V + 1 , k2) 



X- 



Z(S4,N2,k2] 



For large value of N4 and V the previous formula will reduce to the standard 
formula for the average number of baby universe in the canonical ensemble 



A/-rQ4., ^ 5!(N4-V+1)(V+1)Z(S^N4-V,k2)Z(S^V,k2) ,^ ^, 



Z(S4,N2,k2) 



We remember that in each phase the canonical partition function of A four- 
simplices has the behavior ZIS'*, A,k2) x A^'^^'^^^e'^^''^^)^ which substituted in|B]| 
gives 
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A^(S4,N4,k2,V)«^N4(^v(l-^)y^ . (B.7) 

In this way by computing the average number of baby universes with volume 
V in a canonical ensamble of N4 simplices we can obtain informations on the 
susceptibility exponent ystr- This is the reason too for which the argument of the 
baby universes is so popular among people doing numerical simulations. 



Appendix C 
Voronoi Cells 



In this section we shall introduce the definition of Voronoi Cells since we are going 
to use it for the definition of the "metric" dual of a triangulation of a PL-manifold. 

Consider a point set E of a n-dimensional Euclidean space U". We say that 
is discrete if there is a positive real number tq such that Vx,-y G E"^ d(x, y) > 0, 
where d is the distance function between two points of E"^ induced by the Euclidean 
metric [f70[ . 



It follows immediately from the above definition this theorem: 

Theorem C.0.1 IfZ is discrete it has the local finiteness property: for every 
closed ball Br(x) with center m x G E"^ and radius r, Br(x] nZ is a finite set. 

A point set E is relative dense in E"^ if there is a real positive number Rq such that 
all the spheres of E"^ of radius greater than Rq have at least one point of E in their 
interior. Rq is called the covering radius of E. 

A set point E which is discrete and relative dense is a Delaunay or Delone set 



7^ [|71[], and is also indicated as (r, R) from its discrete radius r and covering 
radius R. 

In a Delone set E the r-star, r > 0, of a point x G E is the finite point set Bt-(x) fiE. 

Consider now a Delone set E G E"^ (or, equivalently, a set with a finite number 
of points). Let x G E, the Voronoi cell of x is, by definition, the convex region of 
points of E^ which are more close to x than to any other point of E. 
In practice a Voronoi cell of a point x G E is constructed by considering the 
segments that join x to each point of E in the star of x. Anyway for the moment 
we are considering the case in which E is a sort of periodic crystal so that the 
radius of the star has to be chosen in such a way that it contains x and the points 
of its elementary cell (see figure |C.1|) . 
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X ' > \ 
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Figure C.l: In this regular crystal the points to be considered for constructing the 
Voronoi cell V(x) of x are those ones which belong to the star of x made by the 
points closer to x. These points are also the points which form an elementary cell 
of this crystal 
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Figure C.2: Voronoi cell V(x) of the point x belonging to the crystal of the figure 



Then the n — 1 -dimensional faces of the Voronoi cell of x V(x], called facets, 
are obtained by considering the n — 1 -dimensional orthogonal bisectors to the 
segments. The intersection of the facets V(x) are the n — 2 dimensional face of 
the Voronoi cell and so on (see figure |C.2| ) . 

If we carry on this construction for each point of E we obtain a partition of 
into cells that is called Voronoi tessellation. This tessellation has the following 
properties: 

1) The Voronoi cells are convex region of E"^ and two cells can have in common 
only one n — 1 dimensional facets. 

2) The points of E whose Voronoi cells share a vertex -v lie on a sphere with center 
in "V which has no points of E in its interior. Furthermore it is easy to see that the 
sphere centered in ■y pass among at least n + 1 points of E. 

If E has an infinite number of points and it doesn't have any periodic structure, 
for constructing a Voronoi cell of a point x G E in the general case we need an 
infinite number of operations. But we are going to show that really they are still 
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finite, thanks to the following theorem 

Theorem CO. 2 Let E a Delone set in of type (r, R) in the previous nota- 
tions. To construct a Voronoi cell in of the point x eZ we have to join x 
to all point ofE in the ball B2r(x). 

Dim: 

Consider the Voronoi tessellation of E"^ based on the Delone set E. Let x a point 
of E and V(x) the facets of the Voronoi cell relative to x. By definition the V(x) lie 
on the TL— 1 iperplains that are bisectors of the segment joining x to its neighbor 
points of E. Let 'v any vertex of V(x). We know that ^f is the center of a sphere 
on which there are at least n point of E whose relative Voronoi cells intersect at 
^f. In the interior of this sphere there are no points of E so that d(x,"v) < R. So 
that if we consider the distance between x and any other point y among all other 
at least n points above, we have d(x,-y] < d(x,-v) + d{y,^f] since the reasonings 
for X can be applied equivalently to any other point y, we have d[y,^f) < R and 
then d[x,y] < 2R. This implies that to construct a Voronoi cell in E"^ of a point 
X belonging to a Delone 

set it is sufficient to consider the points of E which are in the ball B2r(x). 

The triangulations I of a PL-manifold, which we are considering, have always 
a finite number of n simplices so that they have a finite number of vertices. On 
I it is defined a metric tensor such that in the interior of each n-dimensional 
simplex it coincides with the Euclidean metric and near the hinges it is like the 
metric tensor of a cone. So it does make sense to define on the simplicial manifold 
the Voronoi cell 

(see reference @ p. 395) of a vertex as the convex region of points of the 



simplicial complex that are closer to the vertex than to any other vertex of the 
simplicial complex. To construct this Voronoi cell relative to a vertex we consider 
on the simplicial complex all the vertices that lie in the star of this vertex and 
the edges joining them with the vertex. The facets of this Voronoi cell lie on 
the n — 1 iperplain orthogonal to the edges in their middle points. Standard 
facts of elementary Euclidean geometry tell us that the vertices of the Voronoi 
cells coincides with the circumcenters of the simplices, that is the centers of the 
spheres circumscribed to the simplices. In this way we have constructed the dual 
of the original triangulation I. In fact to each vertex uniquely it corresponds a 
n-dimensional Voronoi cell, to an edge a n — 1 dimensional facet of the Voronoi 
cell and in general to a k simplex a n — k- Voronoi polyhedron orthogonal to it. In 
particular this dual application maps a n-simplex to its circumcenter. 
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